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(57) Abstract 

This invention is a method and an apparatus of blood measure oxygenation in a subject. A first signal source (106) applies a first 
input signal during a first time interval. A second signal source (108) applies a second input signal during a second time interval. A detector 
(150) detects a first parametric signal responsive to the first input signal passing through a portion of the subject having blood therein ( 102). 
The detector also detects a second parametric signal responsive to the second input signal passing through the portion of the subject. The 
detector generates a detector output signal responsive to the first and second parametric signals. A signal processor (170) receives the 
detector output signal and demodulates the detector output by applying a first demodulation signal to a signal responsive to the detector 
output signal to generate a first output signal responsive to the first parametric signal. Thc'signal processor applies a second demodulation 
signal to the signal responsive to the detector output signal to generate a second output signal responsive to the second parametric signal. 
The first demodulation signal and the second demodulation signal both include at least a first component having a first frequency and a^first 
amplitude, and a second component having a second frequency and a second amplitude. The second frequency is a harmonic of the first 
frequency. The second amplitude is related to the first amplitude to minimize crosstalk from the first parametric signal to the second output 
signal, and to minimize crosstalk from the second parametric signal to the first output signal. 
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METHOD AND APPARATUS FOR DEMODULATING SIGNALS 
IN A PULSE OXIMETRY SYSTEM 

5 Background of the Invention 

Field of the Invention 

The present invention relates to the field of signal processing, and, more particularly, 
relates to the field of processing of signals generated in a physiological monitoring system, such 
as, for example, in a system for measuring blood oxygen saturation using pulse oximetry. 

10 Description of the Related Art 

The present invention will be described herein in connection with a pulse oximetry 
apparatus and a method, which are used to measure blood oxygen saturation in a subject, such as, 
for example, a human patient. The teachings of the present invention can be used in other 
applications wherein useable signal information is obtained in a noisy environment. 

15 In an exemplary pulse oximetry apparatus and a corresponding method, blood oxygen 

saturation is determined by transmitting pulses of electromagnetic energy through a portion of a 
subject which has blood flowing therein (e.g., through a finger, through an ear lobe, or other 
portion of the body where blood flows close to the skin), in the examples described herein, the 
pulses of electromagnetic energy comprise periodic pulses of red light having wavelengths of 

20 approximately 660 nanometers, for example, and periodic pulses of infrared light having 
wavelengths of approximately 905 nanometers. As described, for example, in U.S. Patent No. 
5,482,036 and in U.S. Patent No. 5,490,505, the pulses of red light and the pulses of infrared 
light are applied with the same periodicity but in an alternating and non-overlapping manner. In 
particular, in preferred embodiments, the red pulses are active for approximately 25% of each 

25 cycle and the infrared pulses are also active for approximately 25% of each cycle. The red pulses 
are separated in time from the infrared pulses such that both pulses are inactive for approximately 
25% of each cycle between a red pulse and the next infrared pulse and both pulses are inactive 
for approximately 25% of each cycle between an infrared pulse and the next red pulse. 
(Although described herein below in connection with pulses having 25% duty cycles, it should 

30 be understood by persons of skill in the art that the duty cycles of the pulses can be changed in 
some applications.) After propagating through the portion of the subject, the red pulses and the 
infrared pulses are detected by a detector which is responsive to light at both wavelengths and 
which generates an electrical signal which has a predictable relationship to the intensity of the 
electromagnetic energy incident on the detector. The electrical signal is processed in accordance 

35 with the present invention to provide a representation of the blood oxygen saturation of the 
subject. 
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In conventional time division multiplexing (TDM) demodulation that uses rectangular 
waves to drive the red and infrared LEDs, the conventional process of demodulation using square 
waves can result in the aliasing of the ambient noise components that come close to the 
s.debands of harmonics and the fundamental frequency of the rectangular waves, and the noise 
components are thus collapsed into the output signal generated by the demodulation In 
partrcular, it is very difficult to avoid including hannonics of the line frequency in the 
demodulated output signal. 

Summary of the invention 
The present invention avoids the problems associated with conventional demodulation 
and separation of TDM signals. In particular, the present invention avoids the problem of 
ahasmg of the ambient noise into the passband of the system by selectively demodulating certain 
harmonics of the TDM signal. For example, in one embodiment, only two harmonics (e g the 
fundamental and the first harmonic) are demodulated. In other embodiments, more hannonics 
are demodulated. The present invention specifical.y addresses solutions to problems caused by 
crosstalk resulting from filtering and also resulting from demodulating with only certain 
harmonics instead of demodulating with all harmonics as is done using conventional square wave 
demodulation. I„ a digital implementation of the present invention, the output of the 
photodetector is initially sampled at a very high frequency (e.g., 46.875 Hz), and the signals are 
decmated ( where decimation is lowpass filtering followed by sample rate compression) such 
that the final output signals are generated at a relatively low sampling rate (e.g., 62.5 Hz) which 
provides increased resolution at the output. Thus, bandwidth is traded for resolution in the output 
signal, thus increasing the signal to noise ratio. 

One aspect of the present invention is an apparatus for measuring blood oxygenation in a 
subject. The apparatus comprises a first signal source which applies a first input signal during a 
first tune interval. A second signal source applies a second input signal during a second time 
-nterval. A detector detects a first parametric signal responsive to the first input signal passing 
through a portion of the subject having blood therein. The detector also detects a second 
parametric signal responsive to the second input signal passing through the portion of the subject 
The detector generates a detector output signal responsive to the first and second parametric 
s.gnals. A signal processor receives the detector output signal. The signal processor 
demodulates the detector output signal by applying a first demodulation signal to a signal 
responsive to the detector output signal to generate a first output signal responsive to the first 
parametnc signal and by appiying a second demodulation signal to the signa. responsive to the 
detector output signal to generate a second output signal responsive to the second parametric 
s.gnal. Each of the first demodulation signa. and the second demodulation signal comprises at 
least a first component having a first frequency and a first amplitude and a second component 
havmg a second frequency and a second amplitude. The second frequency is a harmonic of the 
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first frequency. The second amplitude is selected to be related to the first amplitude to minimize 
crosstalk from the first parametric signal to : .he second output signal and to minimize crosstalk 
from the second parametric signal to the first output signal. In one embodiment, the second 
amplitude is determined by turning off one of the first and second signal sources and measuring 
5 the crosstalk between one of the parametric signals and the non-corresponding output signal 
while varying the second amplitude. A second amplitude is selected that minimizes the 
measured crosstalk. 

Another aspect of the present invention is a method of minimizing crosstalk between 
two signals generated by applying a first pulse and a second pulse to measure a parameter. The 

1 o first pulse and the second pulse are applied periodically at a first repetition rate defining a period. 
The first pulse is generated during a first interval in each period, and the second pulse is 
generated during a second interval in each period. The second interval is spaced apart from the 
first interval. The first and second pulses produce first and second parametric signals responsive 
to the parameter. The first and second parametric signals are received by a single detector that 

15 outputs a composite signal responsive to the first and second parametric signals. The method 
comprises the step of applying a first demodulation signal to the composite signal to generate a 
first demodulated output signal wherein the first demodulation signal comprises at least a first 
component having a first frequency corresponding to the first repetition rate. The first 
.component has a first amplitude. The first demodulation signal further comprises a second 

20 component having a second frequency that is a harmonic of the first frequency. The second 
component has a second amplitude which has a selected proportional relationship to the first 
amplitude. The method further includes the step of applying a second demodulation signal to the 
composite signal to generate a second demodulated output signal. The second demodulation 
signal comprises the first component at the first frequency and the first amplitude and further 

25 comprises the second component at the second frequency and the second amplitude. At least one 
■ of the first and second components of the second demodulation signal has a selected phase . 
difference with respect to the corresponding one of the first and second components of the first 
demodulation signal. The method further includes the steps of lowpass filtering the first 
demodulated output signal to generate a first recovered output signal responsive to the first 

30 parametric signal; and lowpass filtering the second demodulated output signal to generate a 
second recovered output signal responsive to the second parametric signal. 

Preferably, the selected phase difference is n. Also preferably, the first pulse and the 
second pulse are generally rectangular pulses having a respective duty cycle. The rectangular 
pulses comprise a plurality of sinusoidal components including a fundamental component 

35 corresponding to the first frequency and a first harmonic component corresponding to the second 
frequency. The fundamental component has a fundamental component amplitude and the first 
harmonic component has a first harmonic component amplitude. The first harmonic component 
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amphtude is related to the fundamental harmonic component amphtude by a first proportionality 
value. The second amplitude of the second component of the first demodulation signal is related 
to the first amplitude of the first component of the first demodulation signal by a second 
proportronality value which is approximately the inverse of the first proportionality value. 

The method in accordance with this aspect of the invention preferably includes the 
further steps of sampling the composite signal when neither the first pulse nor the second pu.se is 
acttve to obtam a sampled signal; and measuring the sampled signal to determine a noise level of 
the parametric signals. 

In a further embodiment according to this aspect of the present invention, the method 
-0 further includes the steps of performing a transform on the composite signal to generate a spectra 
of the composite signal; sampling the spectra at a plurality of frequencies other than at 
predetermined ranges of frequencies around the first frequency and around harmonics of the first 
frequency; determining an average of the magnitudes of the samp.ed plurality of frequencies- and 
companng the average to a selected thresho.d to determine whether the average magnitude 
15 exceeds the selected threshold. 

Another aspect of the present invention is a method of demodulating a composite signal 
generated by applying first and second periodic pulses of electromagnetic energy to a system 
havmg a parameter to be measured and by receiving signals responsive to the electromagnetic 
energy after having passed through the system and being affected by the parameter being 
* measured. The signals are received as a composite signal having components responsive to the 
first and second pulses. The method comprises the step of app.ying a first demodulation signal to 
the composrte signal to generate a first demodulated signal. The first demodulation signal 
compnses a first component having a first frequency corresponding to a repetition frequency of 
the first and second pulses and comprises a second component having a frequency that is a 
s harmonic of the first frequency. The first component has a first amplitude and the second 
component has a second amplitude. The second amplitude has a predetermined re.atronship to 
the first amplitude. The predetermined rebtionship is selected to cause the first demodulated 
s.gnal to have low frequency components responsive only to the first pulse. The method includes 
the fimher step of lowpass filtering the first demodulated signal to generate a first output signal 
> The first output signal varies in response to an effect of the parameter on the e.ectromagnetic 
energy received from the first pulse. 

Preferably, the method in accordance with this aspect of the invention includes the 
further step of applying a second demodulation signal to the composite signal to generate a 
second demodulated signal. The second demodulation signal has first and second components 
■ corresponding to the first and second components of the first demodulation signal. At .east one 
of the first and second components of the second demodu.ation signal has a selected phase 
relatronship with the corresponding one of the first and second components of the first 
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demodulation signal. The method includes the further step of lowpass filtering the second 
demodulated signal to generate a second output signal. The second output signal varies in 
response to an effect of the parameter on the electromagnetic energy received from the second 
pulse. 

Another aspect of the present invention is a pulse oximetry system that comprises a 
modulation signal generator. The modulation signal generator generates a first modulation signal 
that comprises a first pulse that repeats at a first repetition frequency. The first pulse has a duty 
cycle of less than 50%. The modulation signal generator generates a second modulation signal 
comprising a second pulse that also repeats at the first repetition frequency. The second pulse 
has a duty cycle of less than 50%. The second pulse occurs at non-overlapping times with 
respect to the first pulse. Each of the first and second pulses comprises a plurality of components 
wherein a first component has a frequency corresponding to the repetition frequency and wherein 
a second component has a second frequency corresponding to twice the first frequency. The 
second component has an amplitude which has a first predetermined relationship to an amplitude 
of the first component. A first transmitter emits electromagnetic energy at a first wavelength in 
response to the first pulse; and a second transmitter emits electromagnetic energy at a second 
wavelength in response to the second pulse. A detector receives electromagnetic energy at the 
first and second wavelengths after passing through a portion of a subject and generates a detector 
output signal responsive to the received electromagnetic energy. The detector output signal 
includes a signal component responsive to attenuation of the electromagnetic energy at the first 
wavelength and a signal component responsive to attenuation of the electromagnetic energy at 
the second wavelength. A first demodulator multiplies the detector signal by a first 
demodulation signal and generates a first demodulated output signal. The first demodulation 
signal comprises a first component having the first frequency and having a first amplitude. The 
first demodulation signal also comprises a second component having the second frequency and 
. having a second amplitude. The second amplitude has a second predetermined relationship to 
the first amplitude. The second predetermined relationship is approximately inversely 
proportional to the first predetermined relationship. A second demodulator multiplies the 
detector signal by a second demodulation signal and generates a second demodulated output 
signal. The second demodulation signal comprises a first component having the first frequency 
and having the first amplitude. The second demodulation signal further comprises a second 
component having the second frequency and having the second amplitude. At least one 
component of the second demodulation signal has a selected phase relationship with a 
corresponding one component of the first demodulation signal. Preferably, the selected phase 
relationship is a n phase difference. 
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An0ther embodiment incorporates decimation before demodulation. ,„ yet another 
3 mUlti " Chamel dem0dU ' at0r ' ^ ° r With ° Ut ^-duiation decimation is 

demod ,T T ther emb0diment - " 3daPtiVe a,80rithm iS USCd 10 COntr °' the of pre- 

T^ZT:^ ^ ^ adaptive a, g ori th m may 

ZIT ^ ° f 3 l0W?aSS fi,ter " ^ dCCimat0r ^ «" rate 

provided by a signal rate compressor in the decimator. 

Another embodiment of the invention is a method for selecting a sample rate that 
reduces the interference caused by ambient light. 

Brief Description of th» r>r»»,; nr 

draw r PreSem inVenti ° n bC W ° W b COnneCti °" With * e accompanying 

drawing figures in which: Paying 

Figure 1 ilta^es an exa^ block diagran , of , represe „ a , i()n 

saturation in a subject; oxygen 

Figure 2 illustrates exemp.ary waveforms of the current through the LEDs i„ Figure 1 
and the resulting intensities of the red light and the infrared light generated by the LEDs- 

with th F ' 8Ure 3 illUStrateS 2 b ' OCk ° f ^ ° Vera " pro8 — » ^ in -ordance 

with the present invention; 

Figure 4 illustrates a frequency spectra of the fin* modulation signal for « = 0 1 , 

the DC e 7 h he h ° riZOntal rePrCSentS freqU6nCy ^ thC — ^ -P-nts the energy' in' 
the DC and harmonic components of the signal; 

Fi8 "" 5 - «em Pl ^ speemuo „ „. fc , md ^ 

the tadamenal and hann0 „ ic! of conventiona , 6 „ ^ ^ 

Figu-e 6 , he effect „ fc v „ ue of fl on fc sisM| ^ 

S, (0 responsive module p „ lses as the value of B is v»ed wnife , he infiared 
modulation pulses are off; 

Figure 7 illustrates a preferred embodiment of the present invention implemented in a 
digital processing system; 

i„va„,i„„ Fi8Ure 8 i " UStR ' ttS " * a " ed bl0Ck di,gram ° f dem<K " ,tebn ° f *« I"— 

Figu* 0 i„ oslratts . ^ btock diagnm of ^ mo(juiat . on ^ 
* iiven tion , 

Figuae „ musses *. re d dHve wavefo™ ,„ d „ infrared drjve ^ 
by the modulation portion of Figure 9; 
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Figure 1 1 illustrates the demodulation waveforms generated by the demodulation portion 
of Figure 8; 

Figure 12 illustrates a method of time domain sampling the digital detection signal 
during the times when both the red pulses and the infrared pulses are off to obtain information 
5 regarding the level of ambient noise; 

Figure 13 illustrates a block diagram of a system that performs the time domain 
sampling of Figure 12; 

Figure 14 illustrates a method of frequency domain sampling to determine the noise 
floor at frequencies other than the signal frequencies; 
10 Figure 15 illustrates a block diagram of a system that performs the frequency domain 

sampling of Figure 14; 

Figure 16 illustrates a block diagram of the overall processing system in accordance with 
a pre-demodulation decimation embodiment of the present invention; 

Figure 1 7 illustrates a block diagram of a multi-channel processing system in accordance 
1 5 with a pre-demodulation decimation embodiment of the present invention; 

Figure 18 illustrates a block diagram of an adaptive multi-channel processing system in 
accordance with a pre-demodulation decimation embodiment of the present invention; 

Figure 19 illustrates a flowchart of a method for choosing the modulation frequency and 
decimation rate in order to minimize the affects of ambient light; and 
20 Figure 20 is a graph to be used in connection with graphical method for designing a 

demodulation system to minimize interference due to ambient light. 

Detailed Description of the Preferred Embodiment 
Figure 1 illustrates an exemplary block diagram of a representation of a signal 
processing system 100 in accordance with the present invention used to determine blood oxygen 
25 saturation in a subject, such as, for example, a human subject. In the example presented, the 
measurements are performed on a portion of the subject, such as a finger 102 illustrated in Figure 
I. An LED modulation circuit 104 drives a pair of back-to-back light emitting diodes (LEDs) 
106, 108 by applying a periodic signal to the two light emitting diodes 106, 108. The LED 106 is 
selected to emit electromagnetic energy in the red visible light range, and has a wavelength of, 
30 for example, approximately 660 nanometers. The LED 108 is selected to emit electromagnetic 
energy in the infrared range, and has a wavelength of, for example, approximately 905 
nanometers. The LED modulation circuit 104 supplies current in alternating directions so that 
the two LEDs 106, 108 are activated one at a time. In particular, as illustrated by a current 
waveform 120 in Figure 2, current is first applied in a forward direction with respect to the red 
35 LED 106 during a first time interval 122 having a duration x. Thereafter, no current is applied to 
either LED during a second time interval 124 having a like duration x. Then, current is applied in 
a forward direction with respect to the infrared LED 108 during a third time interval 126, also 
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avmg a du ra tion , Then, no current is applied to either LED during a fourth tune interva. ,28 
having a like duration , hereafter, the current is again appHed in the forward direction for the 
red LED 106 during a fifth time interval ,30 which corresponds to the first time interval ,22 It 
can be seen that the overall cyde repeats with a period of duration T equal to 4x. The red LED 
.06 e^ light only when the CUITent „ ^ . r ^ ^ ^ ^ ^ ^ ^ ^ 

LED ,06. Thus, as illustrated by a red intensity waveform 132, the red LED 106 emits light as a 
pulse ,34 during the first time interva, ,22 and as a P u,se ,36 during the fifth time interva, ,30 
and so on. The red pulses repeat with a periodicity equal to T. Similarly, the infrared LED ,08 
emits .nfrared light only when the current is apphed in the fonvard direction with respect to the 
■nfrared LED 108. Thus, as iHustrated by an infrared intensity waveform ,40, the infrared LED 
108 emit, infrared light as a pulse 142 during the third interva, 126. A next infixed pu,se .44 
occurs at an interval T after the infrared pulse 142. Thus, the infrared pulses also repeat with a 
penod,cny equal to T. It can be seen that the red pulses and the infrared pu.ses each have a duty 
cycle of 25o/ 0 , and the red pu.ses and the infrared pulses are separated by inters 

of each period r ( i.e.. the beginning of one pu.se occurs an interva, x after the end of the previous 
pulse). v 

As further iHustrated in Figure the electromagnetic energy pu.ses from the red LED 
.06 and the infrared LED ,08 are app.ied to the finger ,02. A detector ,50 is positioned to 
rece IV e the e.ectromagnetic energy after the energy has passed through a portion of the finger 
102. The detector .50 is se.ected to be responsive to both the red light and the infrared ,ight and 
to generate an output signa, .sponsive to the intensity of the energy received from each source 
An exemp.ary current output signal from the detector .50 is represented by a waveform .52 in 
F.gure 2. As iHustrated, the detector signa, waveform .52 comprises a first pu.se .54 responsive 
to the first red pu.se ,34, a second pu,se ,56 responsive to the infrared P u,se ,42 and a third 
pu* 158 responsive to the second red pu,se ,36. During the time between the first pu,se ,54 
andthesecondpu.se 1 56, the detector signa. waveform .52comprisesno.se . 60, and during the 
tune between the second pu.se .56 and the third pu.se .58, the detector signa. waveform 150 
compnses noise ,62. Tne signa, puises ,54, ,56 and ,58 a,so include noise superimposed 
thereon. Although shown as repeating noise, it shou.d be understood that the noise varies with 
tune. F or example , noise ^ by fflbjen , ^ ^ ^ ^ & ^ ^ 

the 50 Hz or 60 Hz power frequency and their harmonics, P articu,ar,y when the ambient ,ight is 
proved by fluorescent Hghts which generate significant noise at the first harmonic (i.e ,00 Hz 
or ,20 Hz) and the third harmonic (i.e., 200 Hz or 240 Hz). 

output of the detector ,50 is app.ied as an input to a signa, processor b.ock ,70 
which processes the detector signa, and generates a first signa, 5, (,) responsive to the detected 
intensity of the red light incident on the detector ,50 and generates a second signa. *,(,) 
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responsive to the detected intensity of the infrared light incident on the detector 150, As 
illustrated, the signal processing block 170 is synchronized with the LED modulator 104 via a set 
of control lines 180. As will be discussed below, the control lines 180 advantageously 
communicate signals which provide timing information that determines when to activate the red 

, 5 LED 106 and when to activate the infrared LED 108. 

Figure 3 is a pictorial representation of a model of an exemplary system which 
incorporates the present invention. The red LED 106 provides a light intensity represented as 
Ird, and the infrared LED 108 provides a light intensity represented as Iir. The effects of turning 
the LEDs 106, 108 on and off on periodic bases are modeled by a first multiplier or modulator 

10 190 which applies a first modulation signal M\(t) to the red light intensity to generate a 
modulated red signal IrdmodW and by a second multiplier or modulator 192 which applies a 
second modulation signal M 2 (t) to the infrared light intensity to generate a modulated infrared 
signal IirmodW- The modulated light red signal and the modulated infrared signal are applied to 
the finger 102, or other body portion, as described above. The finger 102 has blood flowing 

15 therein and is represented in Figure 3 as a block 102. The blood in the finger 102 has a volume 
and scattering components which vary throughout each cardiac cycle. The blood carries oxygen 
and other materials therein. The oxygen content is a function of both the blood volume and the 
concentration of the oxygen in the blood volume. The concentration of the oxygen in the blood 
volume is generally measured as blood oxygen saturation for reasons which are described in full 

20 in the above- identified issued U.S. Patent Nos. 5,482,036 and 5,490,505. As further described in 
the two referenced patents, the blood oxygen saturation is determined by comparing the relative 
absorption of the red light and the infrared light in the finger 102. The comparison is 
complicated by the noise caused by movement, ambient light, light scattering, and other factors. 

In Figure 3, a pair of signals S { (t) and S 2 (t) represent the effect of the time-varying 

25 volume and scattering components of the blood in the finger 102 on the red light and the infrared 
light, respectively, passing through the finger 102 from the LEDs 106, 108 to the detector 150. 
The red light signal portion S\(t) is caused by the variable attenuation of the red light passing 
through the finger 102. The infrared light signal portion S 2 (t) is caused by the variable 
attenuation of the infrared light passing through the finger 102. To show the effect of the 

30 variable attenuations, the signal portion S\(t) is illustrated as being applied to a first attenuation 
modulator 191 which multiplies the signal by the modulated red output IrdmodW of the first 
modulator 190. Similarly, the infrared light signal portion S 2 (t) is illustrated as being applied to a 
second attenuation modulator 193 which multiplies the signal S 2 (t) by the modulated infrared 
output ItRMODM of the second modulator 192. The outputs of the first and second attenuation 

35 modulators 191, 193 are provided to the receiving photodetector 150. The photodetector 150 is 
modeled as adder 194 and an adder 196. 
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The outputs of the first and second attenuation modulators 1 9 1 , 1 93 are provided to the 
adder 194 to generate a composite signal M(t) where: 

Ht) = S,(0A/,(r) + StfWM. ( i ) 

The signal A/(r) from the adder 194 is provided to the adder 196 where the signal M(t) is 
added to a signal n(t) which represents a composite noise signal caused by ambient light 
electromagnetic pickup, and the like, which are also detected by the photodetector 150 The' 
output of the adder 196 is a signal M(t) = M(t) + „(,) which includes noise components as well as 
the s,gnal components. The noise components include DC components and harmonics of the 
power line frequency that appear in the ambient light. In addition, as will be discussed in more 
detail below, the signal Af«) may also include noise at higher frequencies caused, for example, 
by other devices such as electrocauterization equipment, or the like. 

The M«) signal output of the third adder 196 (i.e., the output of the detector 150) is 
applied to the input of the signal processing block 170. Within the signal processing block 170 
the signal M(t) is first passed through a fixed gain amplifier 197 and then through an analog 
bandpass filter 198. The analog bandpass filter 198 has a passband selected to pass signals in the 
range of 20 Hz to 10,000 Hz. Thus, the analog bandpass filter 198 removes a significant portion 
of the no,se below .0 Hz. The signal components responsive to the blood oxygen saturation are 
frequency shifted by the operation of the two modulation signals A/,(0 and M(0 and are passed 
by the analog bandpass filter 198. 

In the preferred embodiment, the output of the analog bandpass filter 198 is sampled by 
an analog-to-digital converter 199 and converted therein to digital signals. For example the 
s.gnals are preferably sampled at 46,875 samples per second. The output of the ana.og-to-digita. 
converter 1 99 is a signal MF(k). 

The signal MF{k) is provided as a first input to a first demodulating multiplier 210 The 
signal MF(k) is also provided as a first input to a second demodulating multiplier 212 A first 
demodulating signal £>,(*) is provided as a second input to the first demodulating multiplier 210 
and a second demodulating signal Z> 2 (A) is provided as a second input to the second demodulating 
mumpher 212. The output of the first demodulating multiplier 2 1 0 is provided as an input to a 
first lowpass filter 220, and the output of the second demodulating multiplier is provided as an 
mput to a second lowpass filter 222. The bandwidths of the lowpass filters 220, 222 are 
preferably approximately 1 0 Hz. 

The output of the first lowpass filter 220 is a signal S, (/) , which, as discussed below, is 
an estimate of the signal The output of the second lowpass filter 222 is a signal $,(/) 
which, as discussed below, is an estimate of the signal S 2 (t). As will be shown bej the 
select.™ of the first demodulating signal D,(k) and the second demodulating signal D 2 (k) in 
accordance with the present invention substantially reduces or eliminates the effects of noise in 
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the two output signals S^t) and S 2 (t) and also substantially reduces or eliminates crosstalk 
between the two signals. 

In the preferred embodiment of the present invention, the sample rates of the outputs of 
the lowpass filter 220 and the lowpass filter 222 are compressed by respective sample rate 
s compressors 22 1 and 223. In particular, the sample rate compressors 22 1 , 223 reduce the sample 
rate by 750 to a sample rate of, for example, 62.5 Hz to provide an output which can be further 
processed in accordance with the methods and apparatuses described in the above-referenced 
patents. The sample rate compressions which occur in the sample rate compressors 221, 223 

reduce the rate at which the output signals 5,(0 and S 2 (t) need to be processed while 
10 maintaining the sample rate well above the 0-10 Hz frequency content of the signals of interest. 

The outputs of the filters 220, 222, or the sample rate compressors 221, 223, if included, are 

provided on respective output lines 224 and 226. 

In order to facilitate an understanding of how the present invention operates in 

demodulating the output signal MF(k) from the analog-to-digital converter 199, the modulation 
15 signals M,(f) and M 2 U) will first be described in terms of their frequency components. One 

skilled in the art will appreciate that the modulation signals M x (t) and M 2 (t) can each be 

represented as a Fourier cosine series expansion (e.g-,^^ a„cos(rao0, where co = InIT) 

representing the fundamental and harmonic frequencies of the rectangular signal pulses. One 
skilled in the art will understand that the Fourier series expansion includes phases; however, by 
20 suitably selecting the time origin, the phases are set to zero. A component which is 180° out of 
phase with a corresponding component will advantageously be represented by a minus sign 
before the coefficient. 

Figure 4 illustrates a frequency spectra of the first modulation signal M\(t) for n = 0, 1 , 2, 
. . ., where the horizontal axis represents frequency, with the energy in the DC component along 

25 the vertical axis and increasing harmonics of the fundamental frequency along the horizontal 
axis. The length of each component of M\(i) along the vertical axis represents the energy E(n) in 
each component of the frequency spectra. The first component to the right of the vertical axis is 
at the fundamental frequency (i.e., 1/7), which is designated herein as / 0 ; however, it should be 
understood that the fundamental frequency / 0 corresponds to n = 1 . The second component to the 

30 right of the vertical axis is the first harmonic/i (i.e., n = 2), which has a frequency which is twice 
the fundamental frequency. The third component to the right of the vertical axis is the second 
harmonic / 2 (i.e., n = 3), which has a frequency which is three times the fundamental frequency. 
The components to the right of the second harmonic are numbered accordingly. (Note, other 
conventions identify the fundamental frequency as the first harmonic, and designate the second 

35 harmonic as the frequency that is twice the fundamental frequency. The identification of the 
fundamental frequency as^ is used in the discussion that follows.) 
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In Figure 4, a modulation envelope 230 is shown in dashed lines. The modulation 
envelope 230 represents the magnitudes of the fundamental and the harmonics of the signal 
M } (t). The shape of the envelope is determined by the modulation signal Af,(/) which, for a 
repeating rectangular pulse train starting at time / = 0 and having a normalized amplitude of 1, 
can be expressed as: 



M,(t) - I gsinc^)cos(^) (2) 



Where sine is the function (sin 7ix)/7tx (i.e., sinc(«t/7) = sin(wtT/7)/(mtT/7)). In the example 
shown, t = V*T. (Note that for sampled signals, the envelope is more accurately represented as 
sincc/sinP; however, as well known in the art, for the frequencies of interest, the sine function is a 
suitable approximation.) Thus, the frequency spectra has nulls at n = 4, n = 8, n = 12, and so on, 
corresponding to the third harmonic/ 3 , the seventh harmonic / 7 , the eleventh harmonic /„, and so 
on. Note that Equation 2 is an idealized form of the equation for A/,(f), and that in general: 

where a n is a complex number. In the discussion that follows, the values of a„ are assumed to be 
real numbers only. 

A similar frequency spectra (not shown) for the modulation signal M 2 (t) is determined 
by the expression: 



r^f)*Jffl ( 4) 



An envelope for the frequency spectra of second modulation signal M 2 (t) will have the 
same magnitudes; however, it should be understood that because of the (-1)" term in the 
expression for M 2 (t) 9 the fundamental/o and every even harmonic (i.e.,/ 2 ,/ 4 , etc.) are 1 80° out of 
phase with the corresponding harmonic of the first modulation signal M\(t). 

In Figure 3, the analog-to-digital converter 199 converts the signal M (/)to a sequence of 
sampled digital values MF(k) at a sampling rate of, for example, 46,875 samples per second. As 
discussed above, the first demodulating multiplier 210 multiplies the output MF{k) of the 
converter 199 by the first demodulating signal £>,(*) to generate the first output sequence S x (k) , 
and the second demodulating multiplier 212 multiplies the output MF(k) by the second 
demodulating signal D 2 (k) to generate the second output sequence S 2 (k) . The multiplication by 
the multipliers 210, 212 can also be expressed as follows: 

S,(*)«UW*)A<*)] (5) 

and 

S 2 (k) =L?[MF(k)D 2 (k)] (6) 
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where LP is the transfer Function of the lowpass filter 220 and of the lowpass filter 222. 
If, for simplicity, the noise is assumed to be zero, then: 

M(t) = S,(/)A/,(/) + S 2 (/)A/ 2 (/) (7) 

Therefore: 

5 (*) - LP[ [S x (k}Mm + S 2 {k)M 2 {k)]D>(k) ] (8) 

and thus 

S, (k) = LP[ [StikWikWik) + [5 2 (*)M 2 (A)]D,(*) ] (9) 

Similarly: 

S 2 (k) ^L?[[S 2 (k)M 2 (k)]D y (k)^[S l (k)M ] ^D 2 (k)) (10) 
10 Since LP is a linear operator, the right-hand side of Equations 9 and 10 can be split 

into two terms. The first term on the right-hand side of each of Equations 9 and 10 above is 
the desired signal portion of the equation, and the second term on the right-hand side of each 
of the equations is the crosstalk portion. Thus, in order to reduce the crosstalk to zero, the 
second term of each of Equations 9 and 10 is set to zero: 
15 LP[^)M 2 (*)Z),(/:)] = 0 (11) 

and 

L?[S l (k)MmD 2 (k)} = 0 (12) 
By setting the second terms to zero, Equations 9 and 10 reduce to: 

S } (k) =LP[ StkWMDM] (13) 

20 and 

S 2 (k) = LP[ S 2 (k)M 2 (k)D 2 (k) ] (14) 
One goal of the present invention is to select the demodulating signals D } (k) and D 2 (k) to 
satisfy Equations 11 and 12 to thereby reduce Equations 9 and 10 to Equations 13 and 14. This 
is accomplished by utilizing Equations 2 and 3 to simplify the two equations by selectively using 
25 components of the two modulating signals M\(t) and M 2 (t) to generate the demodulating 
sequences D x {k) and D 2 (k). 

In order to simplify the discussion, Equation 2 can be rewritten as: 

OO 

Mi(0 = ££fajcosfwDf> (15) 

where E(>?) is the sine envelope for the fundamental frequency f 0 (n = 1 ) and the harmonics f x (n 
= 2),/ 2 (n = 3), and so on, where cos(mof) represents the cosine term cos(2nnt/T) y where <o=2n/T. 
30 (Note, as discussed above, for discrete sampled signals, the actual envelope of E(«) is a sina/sinp 
function; however, for the frequencies of interest, the sine function is a suitable representation.) 
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As discussed above, the DC term (* = 0) does not need to be considered because of the 
operation of the filter 198, and the analog-to-digital converter 199, as well as the action of the 
demodulation, which shift any unwanted DC or low frequency signals having a frequency less 
than approximately 10 Hz (hereinafter near-DC signals) to higher frequencies before iowpass 
filtering. As a further simplification, the magnitude of the fundamental term in Equation 15 is 
normalized to a value of 1 (i.e., E(l) = 1 ). Note that the normalization results in the need for a 
scale factor, which will be discussed below. Thus, Equation 1 5 becomes: 

M(0 = cos (0/ + a cos 2co/ + b cos 3co/ + c cos 4co/ + . . . (i 6 ) 
The demodulation signal £>,(/) is defined as: 

D t (t) = cos to/ + B cos 2mt ( j ^ 

For reasons set forth below, only the first two cosine terms are needed. 
Similarly, the second modulating signal M 2 (t) becomes: 

M 2 (t) = -cos cor + a cos 2co/ - b cos 3co/ + c cos 4co/ + . . . (jg) 
and the second demodulating signal £> 2 (/) is defined as: 

D 2 (t) = -cos co/ + B cos 2co/ ^ 
Note that the signs of the fundamental and odd harmonics in Equation 18 are 180° out of phase 
with the corresponding terms in Equation 16. 

Note, as will be developed more fully below, by including only the fundamental (cos cot) 
and the first harmonic (cos 2a,/) in each of the demodulation signals, only the signals proximate 
to the fundamental and first harmonic need to be considered. By eliminating higher harmonics, 
the effects of the higher harmonics of the power line frequency are also eliminated in the outpui 
signals generated by the present invention. 

Assume that the filter 198 and the analog-to-digital converter 199 do not affect the 
magnitude of the signal MF(k) with respect to M(t) for the frequencies having significant energy. 
Therefore, starting with Equation 7 above, M(t) can be written as: 
Af(/) = ^(Ofcos tor + a cos 2co/ + b cos 3co/ + . . .] 

+ -StWi-cos cot + a cos2co/ - b cos 3co/ + . . .] ^0) 
When the first demodulating multiplier 210 multiplies A/(/) by D t (t), the terms on the right-hand 
s.de of Equation 20 are multiplied by the terms on the right-hand side of Equation 1 7. Thus: 
M(t)D,(t) = S,(/)[cos co/ + a cos 2tor + b cos 3co/ + . . .][cos co/ + B cos 2co/] 

+ S 2 (/)[-cos cor + a cos 2co/ -b cos 3co/ + . . .][cos <o/ + B cos 2cor] (2 1 ) 

The term S,(0[cos co/ + a cos 2co/ + b cos 3co/ + . . .][cos co/ + B cos 2co/] is the signal term which 
is to be preserved, and the term 5 2 (/)[-co s co/ + a cos 2co/ -b cos 3co/ + . . .][cos co/ + B cos 2m] is 
the crosstalk term to be eliminated. 
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Expanding the crosstalk term from Equation 21, generates: 

crosstalk = S 2 (/)[-cos 2 co/-B cos co/ cos 2co/ + a cos2co/ cos co/ + aB cos 2 2co/ 

- b cos 3co/ cos cot - bB cos 3co/ cos 2co/ + . . .] (22) 
Using the identity, cos(x)cos(y) = V4[cos(x+y)+cos(x-y)], the crosstalk term from 

- 5 Equation 22 becomes: 

crosstalk = S 2 (t)[ -!4(cos 2co/ + 1 ) + ((a - B)/2)[cos 3co/ + cos co/] 
+ (aB/2)[cos 4co/ + 1 ] - (b/2)[cos 4co/ + cos 2co/] 

- (bB/2)[cos 5co/ + cos co/] + . . .] (23) 
The remaining terms in Equation 23 will all have a factor of cos co/ or higher. Thus, Equation 23, 

10 when fully expanded only includes near- DC terms: 

crosstalk D c = LP[ 1 S , 2 (/)[(aB/2)- , / 2 ]] (24) 
where S 2 (t) corresponds to the infrared portion of the original plethysmograph signal which has a 
bandwidth of interest of approximately 0 to 10 Hz. Any components present above 1 0 Hz will be 
eliminated by the action of the lowpass filter 220. Thus, it can be seen that only the signals of 

15 interest are folded back to DC or near-DC. By using the lowpass filter 220, the DC terms and 
near-DC terms can be isolated so that only the DC terms and near-DC terms of the crosstalk are 
presented at the output of the lowpass filter 220. Thus, in order to eliminate the crosstalk, the 
crosstalk terms in Equation 24 need to be set to zero: 

LP[S 2 (/)[aB/2-/ 2 ]] = 0 (25) 

,20. Thus: 

B=l/a (26) 
The result in Equation 26 can also be expressed using a geometric interpretation of 
vector projection (i.e., dot products) of &(/) and S,(/) wherein the projection of S 2 (t) onto D,(/) is 
equal to zero and the projection of S 2 (t) onto D 2 (t) is maximized. In other words, express S,(/), 

25 S 2 (t), AO) and D 2 (t) as vectors of samples in an w-dimensional sample space (e.g., S,(/) is 
represented as a vector Si of samples S\(k)). For example, in a preferred embodiment, n = 148, 
and thus Si, S 2 , D, and D 2 are vectors of 148 samples each. The first crosstalk term is S r D 2 . The 
second crosstalk term is S 2 Di. The first signal output is S r D,. The second signal output is 
SvD 2 . Select the vectors Di and D 2 to drive the crosstalk terms to zero. 

30 The relationship in Equation 26 also works to preserve the signal term. In particular, the 

signal term in Equation 21 can be expanded and lowpass filtered in the same manner as the 
crosstalk term to obtain: 

signal = S t (/) = LP[S,(/)[(aB/2) + 14]] (27) 
Using the relationship from Equation 26, then Equation 27 becomes: 

35 signal = S t (0 = LP[S,(/)[(a/2a) + 14] = LP[5,(/)] = S,(i) (28) 
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It can be readily shown that the same relationship holds for the crosstalk term and the 
signal term for the signal S 2 (t) by defining the second demodulation signal D 2 (t) as: 

A>(0 = -cosa>r+Bcos2<»/ (29) 
and multiplying M 2 (f) by D 2 (t). After expanding the crosstalk and signal terms and eliminating 
the terms above 10 Hz, it can be shown that by selecting B = 1/a. the crosstalk term is canceled 
and the signal term $(/) is recovered. 

From the foregoing, it can be seen that by choosing the relationship between the 
magnitude of B as the reciprocal of a, then the crosstalk terms are eliminated and the signal lerms 
are preserved. Note that neither A nor B is an absolute value. As set forth in Equation 16, a is 
the magnitude of the cos?©* term of M,(/) when the magnitude of the cos©/ term of A/,(r) is 
normalized to 1. Similarly, from Equation 17, B is the magnitude of the cos 2at term of D t (t) 
when the cos cor term of D\{t) is normalized to 1 . 

It should be understood that both D,(0 and D 2 (t) can include higher harmonic terms; 
however, such additional terms could result in increased sensitivity to the noise of fluorescent 
lights and the like because of the harmonics of the 60 Hz power line frequency (or the 50 Hz 
power line frequency in other countries). For example, Figure 5 illustrates an exemplary 
spectrum of the first and second harmonics of the present invention when the fundamental 
frequency is selected to be 3 16.7 Hz. Thus, the first harmonic frequency is 633.4 Hz. Note that 
the variations in the signals caused by blood flow throughout a cardiac cycle causes the 
fundamental and harmonics of modulation frequency to be surrounded by sidebands representing 
the frequency content of the plethysmograph. For example, in Figure 5. the first and second 
harmonics are at 3 16.7 Hz and 633.4 Hz, ±10 Hz. 

As further illustrated in Figure 5, the conventional 60 Hz power line frequency has 
harmonics at 120 Hz, 180 Hz, 240, etc. Thus, the nearest harmonics of the power line frequency 
to the first harmonic of the present invention are at 300 Hz and 360 Hz, and the nearest 
harmonics of the power tine frequency to the second harmonic of the present invention are 600 
Hz and 660 Hz. Similarly, if used in a country having a 50 Hz power line frequency, the nearest 
harmonics to the first harmonic of the present invention are 300 Hz and 350 Hz, and the nearest 
harmonics to the second harmonic of the present invention are 600 Hz and 650 Hz. Even if the 
power frequency were to vary by up to 1 .5 percent, the noise generated by the ambient light from 
fluorescent lamps, or the like, would not be at the first and second harmonic frequencies of the 
present invention. The fundamental frequency has thus been selected to avoid power line caused 
ambient noise at the first and second harmonic frequencies. 

The foregoing discussion assumed that the filter 198 did not significantly affect the 
amplitude of the filtered signal. If the filter 198 does have an affect on the amplitude, then B will 
be a constant times the value of B determined above: 

B = k/a (30) 
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where k depends on the relative attenuation of the first harmonic and the second harmonic 
through the filter 198. 

Although the value of the coefficient B can be calculated as set forth above, the 
calculations may be complicated if the filter 198 or the modulators 190, 192 introduce phase 
5 ^ changes which cause the calculations to be performed on complex numbers. For example, if the 
modulation signals M\ (/) and M 1 (t) are not rectangular waves which have 25% duty cycles and 
which are precisely 1 80° out of phase, as illustrated herein, then the coefficients of the frequency 
components of the modulation signals may be complex to account for the phase relationships, 
and thus, the coefficients of the demodulation signals may be complex. 
10 . As illustrated in Figure 6, the value of B can also be determined empirically by 

performing a initial measurement with one channel (i.e., either the red pulse or the infrared pulse 
turned off) and minimizing the crosstalk. In particular, during the initial measurement, the 
waveform 140 in Figure 2 is set to a continuous zero value so that no infrared pulses are 
generated. Thus, the detector 150 (Figure 1) receives only the light generated by the red LED 

15 1 06. Thus, M 2 (t) is set to zero, and Equation 1 0 for S 2 (t) becomes: 

5 2 (/) = LP[5 1 (/)M,(/)A(0] (31) 
it can be seen that S 2 (/) includes only a crosstalk portion, which can be measured on the output 
from the second lowpass filter 222. Thus, by varying the value B while monitoring the 
.magnitude or the RMS (root-mean-squared) value of the output signal S 2 (t) , a minimum 

20 magnitude S 2 (t) min for the output signal S 2 (t) can be found which corresponds to the best 
value Bbest for B. In an ideal system, the best value for B corresponds to a zero value for the 
output signal S 2 (t) ; however, in a real environment, the best value of B may correspond to a 

non-zero value for S 2 (t) (i.e., a minimum error for S 2 (t) ). It should be understood that the 
value of Bbest can also be determined by turning off the red LED 106 and varying B while 

25 monitoring S, (/) until S } (t) is minimized. 

From the foregoing, it can be seen that the effect of the modulation signals D } (t) and 
D 2 (t) is to shift the DC or near-DC noise terms up in frequency while shifting the signals of 
interest at the harmonics back to DC or near-DC, which in effect interchanges the noise spectra 
and the signal spectra so that the noise spectra can be eliminated by the action of the lowpass 

30 filters 220, 222, leaving only the signals of interest. 

Figure 7 illustrates a preferred embodiment of the present invention which implements 
the functions described above in a digital system. Preferably, the digital system comprises a 
digital signal processor (not shown), and the blocks described herein comprise data structures 
within the digital signal processor and software routines that implement the processes described 
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below. In particular, the present invention comprises an LED demodulation block 300 which 
receives a digital configuration signal on a bus 310, a clock signal on a line 312 and a digital 
detector signal on a bus 3 14 as inputs. The digital configuration signal bus 3 10 provides a way to 
change the configuration of the LED demodulation block 300 to accommodate different LEDs 
and different detection algorithms. Preferably, the clock signal on the line 3 12 is a 46,875 Hz 
(46.875 kHz) square wave signal which is used to synchronize the timing functions of the present 
invention. The digital detector signal on the line 314 is the output of the analog-to-digital 
converter 199. The analog-to-digital converter 199 is connected to the output of the detector 150 
(via the amplifier 197 and the filter 198) and samples the output of the detector 150 at 46,875 
samples per second to provide a stream of sampled digital values of the red light and infrared 
light incident on the detector 150. 

The LED modulation block 300 generates a demodulated red signal output on a bus 340 
and generates a demodulated infrared signal output on a bus 342. The demodulated red signal 
output is passed through the low pass filter 220 and is output therefrom as the signal S, (/) . The 
demodulated infrared signal output is passed through the low pass filter 222 and is output 
therefrom as the signal S 2 (t) . As further illustrated in Figure 8, the LED demodulation block 
300 comprises a modulo-M block 350, an LED demodulation state table block 352, the first 
demodulating multiplier 210 and the second demodulating multiplier 212. 

The modulo-M block 350 receives the main 46,875 Hz clock signal on the line 312 as 
one input and receives a MODULUS signal on a bus 354 as a second input. The bus 354 forms a 
portion of the configuration bus 310. The modulo-M block 350 divides the clock signal by the 
MODULUS signal and generates a RESIDUE signal (described below) on a bus 356 which is 
provided as one input to the LED modulation state table block 352. The LED modulation state 
table block 352 also receives the configuration signals on the configuration bus 3 10. 

The LED demodulation state table is responsive to the residue signal and the 
configuration signals to generate the first demodulating signal Z>,(0 on a bus 360 and to generate 
the second demodulating signal D 2 (t) on a bus 362. The first demodulating signal />,(/) is 
provided as one input to the first demodulating multiplier 210, as described above. The second 
demodulating signal £> 2 (/) is provided as one input to the second demodulating multiplier 212, as 
described above. The first demodulating multiplier 210 and the second demodulating multiplier 
212 receive the digital detector signal on the line 314 as respective second inputs. The 
demodulating multipliers 210, 212 multiply the digital detector signal by the first demodulating 
signal A(r) and the second demodulating signal D 2 (t), respectively, to generate a demodulated 
red signal and a demodulated infrared signal on the buses 340 and 342, respectively. Because the 
outputs of the two demodulating multipliers 210 and 212 include the terms cos or, cos 2cor, and 
higher, the demodulated signals on the buses 340 and 342 are provided as respective inputs to the 
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low pass filters 220 and 222 to pass only the near-DC terms, as discussed above. The outputs of 
the iowpass filters 220 and 222 on the buses 344 and 346, respectively, are the 5, (r) signal and 
the S 2 (t) signal which contain only the near-DC terms, which, in accordance with the 
discussion presented above represent the original input signals S,(r) and S 2 (t) with the unwanted 
5 noise substantially reduced or eliminated. The two signals S,(r) and S 2 (t) are then applied to 
computation circuitry (not shown) which computes the blood oxygen saturation and other 
cardiography parameters in a manner described in the above-cited U.S. Patent Nos. 5,482,036 
and 5,490,505. 

The residue signal generated as the output from the modulo-M block 350 is a multiple 

10 bit signal that counts from 0 to MODULUS- 1. In the preferred embodiment described herein, 
MODULUS has a value of 148. Thus, the RESIDUE output of the modulo-M block 350 counts 
from 0 to 147. The RESIDUE output of the modulo-M block 350 is a number that is provided as 
the input to the LED demodulation state table block 352. As illustrated in Figure 9, the 
RESIDUE output on the bus 356 corresponds to the signal 1 80 in Figure 1 and is also provided to 

15 the input of an LED modulation state table block 370 which, together with an LED driver circuit 
372, comprise the modulation block 104 (Figure 1) which generates the drive signals to the red 
LED 106 and the infrared LED 108. As described above, the red LED 106 and the infrared LED 
108 generate the modulation signals M,(/) and A/ 2 (/), respectively, which effectively operate as 
carriers for the plethysmograph waveform to be measured. In particular, as illustrated by a red 

20 drive timing waveform 374 and by a infrared drive timing waveform 376 in Figure 10, the 
modulation state table block 370 generates a red signal pulse 378 during the time when the 
RESIDUE signal increments from 0 to 36. Then, the modulation state table block 370 generates 
neither a red signal pulse nor an infrared signal pulse during the time when the RESIDUE signal 
increments from 37 to 73. Then, the modulation state table block 370 generates the infrared 

25 signal pulse 380 during the time when the RESIDUE signal increments from 74 to 1 10. Then, 
the modulation state table block 370 again generates neither a red signal pulse nor an infrared 
signal pulse during the time when the RESIDUE signal increments from 111 to 147. The 
RESIDUE signal then resets to 0 and the process repeats continuously. 

The red signal pulse 378 and the infrared signal pulse 380 from the modulation state 

30 table block 370 are provided as inputs to the LED driver circuit 372 which turns on the red LED 
106 when the red signal pulse 376 is active and turns on the infrared LED 108 when the infrared 
signal pulse 378 is active by generating the current waveform 120 illustrated in Figure 2. The 
circuitry for converting the red signal pulse 376 and the infrared signal pulse 378 to the bi- 
directional current pulses of the waveform 120 is conventional and does not need to be described 

35 herein. 
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In the preferred embodiment the LED demodulation state table block 352 implements 
demodulation equations which generally correspond to the Equations 17 and 19 described above, 
in particular, the LED demodulation state table block 352 receives the RESIDUE as one input to 
the state table and steps through the state table based upon the current value of the RESIDUE. 
The LED demodulation state table block 352 generates two output values for each value of the 
RESIDUE, wherein the first output value is the first demodulation signal Z),(0 on the signal bus 
360, and the second output value is the second demodulation signal D 2 (t) on the signal bus 362. 

In particular, the LED demodulation state table block 352 implements the following 
forms of the demodulation signal D\{t) and the D 2 (t) equations: 



D,(t) = -SCL\cos 



2nt 



+HWD 



cos 



V 



4nt 



R - 18 J -HWA\ 
Modulus J 



R - 18.5 -HWA\ 
Modulus J 



V) 



(32) 



and 



D 2 (t) = -SCL 



-HWD 



cos 4nt 



2 (R-18.5-HWAX 
V Modulus J 



18.5 -HWA\ 
Modulus J / 



(33) 



10 



In Equations 32 and 33, the value SCL is a scale factor which determines the magnitudes of the 
two demodulation signals and which is used to compensate for the normalization discussed 
above and to compensate for other factors, such as, for example, non-ideal rectangular pulses. 
The method of determining the scale factor will be set forth below. In one particularly preferred 
embodiment, the value of SCL is 2.221441469. The value HWD is a hardware distortion factor, 
which corresponds to the value of B discussed above. The determination of the value B was 
described above, and will be described again below in connection with this preferred 
embodiment. In one particularly preferred embodiment where the pulses applied to the red LED 
1 06 and the infrared LED 1 08 are idealized rectangular waves having 25% duty cycles, the value 
of HWD can be calculated to be 1.414213562. This ideal value for HWD can be determined by 
recognizing that the value of the coefficient A for the cos 2cof terms in Equations 16 and 18 is 
determined by the sine function. When the coefficient of the cos co/ term is normalized to 1, as in 
the two equations, then the value of the coefficient a is equal to . Thus, the ideal value for 
B (i.e., HWD) is ^2. Of course, the actual value of the coefficient B, and thus HWD, will vary 
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when the red pulses and the infrared pulses are not true rectangular waves. Since, in actual 
embodiments, the pulses will have finite rise times and fall times, the optimum value of HWD is 
preferably found empirically in the manner described below. 

The value 1 8.5 in Equations 32 and 33 is used to align the demodulation waveforms with 
5 , the modulation waveforms so that the peak of the cosine functions corresponds to the midpoints 
of each of the modulation waveforms. The value HWA is a hardware delay factor which may be 
needed in certain embodiments to compensate for delays in the analog processing, the digital 
processing or both, which cause the demodulation signals D\(t) and D 2 (t) to be out of phase with 
the modulation signals M\(t) and M 2 {t). In an ideal environment, the value of the hardware delay 
10 factor is 0. However, in one particularly preferred embodiment, the value of the hardware delay 
factor is 39. The modulus was described above and is basically the number of steps in each 
period of the waveforms. In the embodiment described herein, the modulus is 148. The value R 
is the RESIDUE, which varies from 0 to modulus- 1, and thus, in the preferred embodiment, R 
varies from 0 to 147. 

15 In operation, the clock signal on the line 3 12 causes the moduio-M block 350 to generate 

the RESIDUE signal, as described above. The RESIDUE value is applied to the LED 
modulation block 104 which generates the modulation signals and M 2 (t), as described 
above. The RESIDUE value is also applied to the LED demodulation state table block 352 
which generates a new value for D { (t) and a new value for D 2 (t) for each new RESIDUE value. 

20 Thus, 148 values of /),(/) and D 2 (t) are generated for each complete cycle. Because the clock 
signal is operating at 46,875 Hz, the modulation signals A/,(f) and M 2 {i) and the demodulation 
signals Z),(/) and D 2 (() have a fundamental frequency of 316.722973 Hz, which, as discussed 
above, does not correspond to any harmonic of conventional 50 Hz or 60 Hz power line 
frequencies. 

25 The HWA (hardware delay factor) value, the HWD (hardware distortion factor) value 

and the SCL (scaling factor) value are found empirically as follows. First, the ideal values of the 
hardware delay factor, the hardware distortion factor and the scale factor are applied to the 
Equations 32 and 33 in the LED demodulation state table block 352 (i.e., HWA= 0, HWD = 
1 .414213562, and SCL = 2.221441469). To determine the optimum value of the hardware delay 

30 factor, the second modulation signal M 2 (t) is set to a constant value of zero (i.e., the infrared LED 
is maintained in its OFF state). The red LED pulses are applied as set forth above, and the digital 
detector output signal from the analog-to-digitai converter is monitored and compared to the 
modulation signal M,(/). The relative delay between the beginning of the modulation signal 
M,(/) and the detection of the beginning of the responsive output from the analog-to-digital 

35 converter is the optimum hardware delay factor (HWA) value. In one exemplary embodiment, 
the optimum value of the hardware delay factor is 39. 
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After determining the value of the hardware delay factor and applying it to Equations 32 
and 33, the ideal value of the hardware distortion factor and the ideal value of the scale factor are 
applied to the two equations. Again, with the red LED pulses applied to the red LED 106 and no 
pulses applied to the infrared LED, the value of the hardware distortion factor is slowly varied 
from its ideal vaiue while the DC component of the demodulated infrared signal output on the 
line 342 is monitored. The value of the hardware distortion factor is varied until the measured 
DC component is minimized, and the value of the hardware distortion factor corresponding to the 
minimal DC component is selected as the optimum value for the hardware distortion factor. 

Next, with the value of the hardware delay factor and the value of the hardware 
distortion factor set to their respective optimum values, as determined above, the value of the 
scale factor (SCL) is initially set to t. Again, with the modulation system generating pulses only 
to the red LED 106, the DC component of the demodulated red signal output on the line 340 is 
measured. In addition, the difference in amplitude between the on state and the off state of the 
digital detector signal from the filter 198 is measured. The ratio of the measured amplitude 
difference to the measured DC component of the demodulated red signal output is selected as the 
optimum value for the scale factor. 

An exemplary demodulation waveform D,(t) is illustrated by a waveform 400 in Figure 
1 1 and an exemplary demodulation waveform D 2 (t) is illustrated by a waveform 402 in Figure 
1 1. The demodulation waveforms in Figure 1 1 are illustrated with the hardware delay factor set 
to 0 in order to align the waveforms with the modulation waveforms in Figure 10. It should be 
understood that when the hardware delay factor is non-zero, the demodulation waveforms in 
Figure 1 1 will be shifted in phase with respect to the modulation waveforms in Figure 10. 

Although described above in connection with the variation of the amplitude of the first 
harmonic component of the demodulation signals in order to minimize the crosstalk, it should be 
understood that the relative amplitude of the second harmonic component of the demodulation 
signals with respect to the amplitude of the fundamental component of the demodulation signals 
is determined by the relationship of the amplitude of the first hannonic component of the 
modulation signals to the amplitude of the fundamental component of the modulation signals. 
The relationship of the amplitude of the first hannonic component of the modulation signals 
depends in part upon the duty cycles of the modulation signals. If the modulation duty cycles are 
varied, the amplitude of the first hamionic component of the modulation signals changes. Thus, 
the crosstalk may also be minimized by holding the amplitudes of the components of the 
demodulation signals constant while varying the duty cycles of the modulation signals. One 
skilled in the art will appreciate that other variations in the modulation and demodulation signals 
35 may also be used to minimize the crosstalk between the two output signals. 

A plurality of signals S h S 2 , S, . . . S n can be demodulated and the crosstalk between 
signals reduced to a minimum by application of the foregoing invention to more than two signals. 
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Additional information can advantageously be derived from the digitized detection 
signal on the bus 314 and can be used to provide indications regarding the reliability of the 
demodulated signals generated as described above. In particular, although the present system is 

capable of demodulating the £,(/) signal and the S 2 (t) signal in the presence of significant 
5 ,Br ambient noise from light and other sources, it is possible that the level of the ambient noise is 
sufficiently high to affect the demodulated signals. Figures 12 and 13 illustrate a time domain 
method and system for determining the ambient noise level, and Figures 14 and 15 illustrate a 
frequency domain method and system for determining the ambient noise level. 

As illustrated in Figure 12, the digital detection signal 152 is sampled by a sample signal 

10 represented by a waveform 500, which comprises a plurality of sampling pulses 502. The 
sampling pulses 502 are timed to occur during the intervals between the red pulses 134, 136 and 
the infrared pulses 142, 144 when no red light and no infrared light should be detected by the 
detector 150 (Figure 1). Thus, any energy detected during the sample intervals is primarily 
caused by ambient light and other noise sources. As illustrated, the sampling pulses 502 

15 preferably occur at the approximate midpoint of each interval between the red and infrared 
pulses. 

As illustrated in Figure 13, the digital detection signal bus 314 is provided as an input to 
a time domain sampler 520. The time domain sampler 520 also receives the RESIDUE signal on 
the bus 356 as a second input The time domain sampler is responsive to the RESIDUE signal to 

20 sample the digital detection signal at times when the value of the RESIDUE signal corresponds 
to the quiescent times of the red pulses 134, 136 and the infrared pulses 142, 144. As described 
above, the red pulses 134, 136 are generated when the RESIDUE signal has values between 0 
and 36, and the infrared pulses are generated when the RESIDUE signal has values between 74 
and 1 10. Thus, assuming no hardware delay, the sampling pulses 502 are preferably generated, 

25 for example, when the RESIDUE signal has a value of 55 and when the RESIDUE signal has a 
value of 129, which positions the sampling pulses at the approximate midpoints of the quiescent 
intervals between the pulses. As discussed above, the actual system has a hardware delay caused 
by processing times. Thus, if the system has a hardware delay factor of, for example, 39, the 
sampling pulses 502 are shifted in time to occur when the RESIDUE signal has a value of 94 and 

30 a value of 20 (168 modulo M8 ). The sample times used by the time domain sampler 520 are 
advantageously determined by configuration signals received via the digital configuration bus 
310, described above. For example, the time domain sampler 520 can be initially set to sample at 
RESIDUE signal values of 55 and 129, and the value of the hardware delay value factor (HWA) 
communicated by the digital configuration bus 310 is added to both values to shift the sample to 

35 the correct sample interval. 

As illustrated in Figure 14, a detection signal spectra 550 includes the two frequency 
components corresponding to the fundamental and first harmonic of the modulation signal at 
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316.7 Hz and 633.4 Hz, respectively. The spectra 550 further includes the fundamental and 
multiple harmonics of the 60 Hz power line frequency. In addition, the spectra 550 includes 
noise at a multitude of frequencies which may be caused by various sources. One particularly 
troublesome source of noise encountered in pulse oximetry systems is an electrocauterization 
device, which uses a high frequency electrical current to make surgical incisions and to cauterize 
the surrounding blood vessels at the same time. Although primarily high frequency noise 
sources, such devices also generate significant noise at lower frequencies because of arcing. 
When an electrocauterization device is operated close to a pulse oximeter detector, the noise 
generated by the device can overwhelm the signals generated by the pulse oximetry detector. In 
other words, the noise floor can be greater than the detectable signal from the pulse oximetry 
detector. 

It is desirable to detect when the noise floor is too high so that the pulse oximetry system 
can indicate that the demodulated signals may not be reliable. In order to determine the level of 
the noise floor, the present invention samples the spectra 550 to determine the content of the 
frequency components detected at frequencies other than the fundamental and harmonic 
frequencies of the modulation signals. In particular, as illustrated by a sample control signal 560 
in Figure 1 4, the portions of the spectra 550 which do not include the fundamental and harmonics 
of the modulation signal are sampled. Thus, in the preferred embodiment, the magnitudes of the 
spectra at 316.7 Hz, 633.4 Hz, 950.1 Hz, etc., are not sampled. Furthermore, because a band of 
frequencies around the fundamental and harmonics of the modulation signal also include 
significant information caused by the modulation of the red pulses and the infrared pulses by the 
changes in blood flow during each cardiac cycle. Thus, as illustrated in Figure 14, in the 
preferred embodiment, a band of frequencies surrounding the fundamental and harmonic 
frequencies of the modulation signals (i.e., the sidebands discussed above) are not included in the 
samples. For example, a band of at least ±10 Hz around each of the fundamental and harmonic 
frequencies is not included in the samples. 

The intensities at the sampled frequencies are averaged, and an output signal is generated 
which represents the average intensity of the noise signals. Other portions (not shown) of the 
digital processing system advantageously monitor the average intensity of the noise signals, and, 
if the average intensity exceeds a selected threshold based upon the size of the measured 
plethysmograph, then the demodulated output signals from the system are considered as being 
unreliable and should not be used. 

Figure 1 5 illustrates a prefeired embodiment of a system that determines the noise floor, 
as described above. The system of Figure 1 5 includes a Fast Fourier Transform block 600 which 
receives a plurality of samples from the digitized detector bus 314 and generates a transformed 
output on a bus 610. The transformed output on the bus 610 represents the spectra of the 
samples. In the preferred embodiment, a sufficient number of samples are taken to represent 
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approximately 44 milliseconds of data so that at least two cycles of the 60 Hz power are included 
within the samples. For example, approximately 1,024 samples can be taken during the 44- 
millisecond interval at a sample rate of approximately 23.4 kHz (e.g., one-half the system timing 
rate). The spectra for a 44-millisecond interval are provided as inputs to a spectral sampler 620 
5 .which eliminates the samples in the ±10 Hz bands around the fundamental and harmonic 
frequencies of the modulation signals. The output of the spectral sampler 620 is provided on a 

- bus 630 and is thereby provided as an input to an averager 640. The averager 640 averages the 
sampled noise spectra which it receives and provides an averaged output on a bus 650. The 
averaged output on the bus 650 represents the noise floor and is provided to other portions of the 

10 digital processing system where it is compared to the selected threshold to determine whether the 
noise floor is excessive. The threshold is not necessarily fixed, but is dependent on the strength 
of the plethysmograph, which in turn depends upon the perfusion of blood in the body portion 
being measured. 

The embodiment of Figure 15 can also advantageously be used to determine whether the 
15 ambient noise is primarily at 60 Hz, corresponding to power line frequencies in the United States 
and Canada, or at 50 Hz, corresponding to power line frequencies in Europe. The foregoing 
modulation frequency of 316.7 Hz is selected to avoid the harmonics of the 60 Hz power line 
frequency as well as the 50 Hz power line frequency. If a significant shift in the power line 
frequency is detected such that aliasing of the ambient noise occurs at the frequencies of interest, 
20 then the modulation frequency can be changed to displace the modulation harmonics farther from 
the harmonics of the power line frequency, such as, for example, by changing the 46,875 Hz 
sampling frequency, or by changing the modulus. 
Pre-Demodulation Decimation 

For convenience, the previous embodiments do not show the signal MF(k) being 
25 decimated before demodulation. However, as discussed in more detail below, the signal MF{k) 
can advantageously be decimated prior to demodulation. The pre-demoduiation decimation 

- technique can reduce the computational burden required to perform the demodulation operations, 
primarily because the decimated sample rate is lower than the original (undecimated) sample 
rate. Computation can also be reduced because, as will be seen, the numerical sequences used in 

30 the demodulator are, in some circumstances, shorter than the sequences given in Equations 32 
and 33. Pre-demodulation decimation is a generalization of the previous embodiments and 
reduces to the previous embodiments when the pre-demodulation decimation rate is one. 

Figure 16 is a pictorial representation of a system that incorporates pre-demodulation 
filtering and decimation. Figure 16 is similar to Figure 3, and like numbers refer to like elements 

35 in the two figures. Figure 16 shows the first modulator 191 having a signal input and a 
modulation input M } (t). The second modulator 193 has a signal input S 2 (t) and a modulation 
input M 2 (t). The pair of signals S x (t) and S 2 (t) represent the effect of the time-varying volume 
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and scattering components of the blood in a finger (or other body part) on the red light and the 
infrared light, respectively, passing through the finger. The red light signal portion S t (t) is caused 
by the variable attenuation of the red light passing through the finger 102 (shown in Figure 1). 
The infrared light signal portion S 2 (t) is caused by the variable attenuation of the infrared light 
passing through the finger 102. The outputs of the first and second modulators 191, 193 are 
provided to the receiving photodetector 150. The photodetector 150 is modeled as an adder 194 
and an adder 196. The outputs of the first and second modulators 191, 193 are provided to the 
adder 194 to generate a composite signal M(t) where: 

MO = S,(/)M ( (0 + S 2 (/W 2 (0. (34) 
The output signal M(t) from the adder 194 is provided to an adder 196 where a signal 
«(0 is added to the signal M(t). The signal »(/) represents a composite noise signal caused by 
ambient light (including DC and harmonics of the power line frequency), electromagnetic 
pickup, and the like, which are also detected by the photodetector 1 50. In addition, the signal n( t ) 
may also include noise at higher frequencies caused, for example, by other devices such as 
electrocauterization equipment, or the like. The output of the adder 196 is a signal Af(t) = M{t) + 
"(0 which includes noise components as well as the signal components. 

The M(t) signal output of the adder 196 (i.e., the output of the detector 1 50) is applied to 
the input of a signal processing block 1600. Within the signal processing block 1600, the signal 
M(t) is first passed through the amplifier 197 and then through the analog bandpass filter 198. 
The analog bandpass filter 198 provides anti-aliasing and removal of low frequency noise and 
DC. The filter 198 has a passband selected to pass signals in the preferred range of 20 Hz to 
10,000 Hz. The analog bandpass filter 198 removes a significant portion of the noise below 20 
Hz. The signal components responsive to the blood oxygen saturation are frequency shifted by 
the operation of the two modulation signals M,(/) and M 2 (t) and are passed by the analog 
bandpass filter 198. 

In one embodiment, the output of the analog bandpass filter 198 is sampled by the 
analog-to-digital converter 199 and converted therein to digital signals. In one embodiment, the 
signals are sampled at 46,875 samples per second. The digital signals from the analog-to-digital 
converter 199 are provided as inputs to a lowpass digital filter 1620. Output signals from the 
digital filter 1620 are provided to a sample rate compression block 1622 that reduces 
(compresses) the sample rate by a decimation rate The lowpass digital filter 1620 and sample 
rate compressor 1622 together comprise a decimator 1621 (decimation comprises lowpass 
filtering followed by sample rate compression). The digital filter 1620 provides anti-aliasing 
filtering and the sample rate compression block 1622 preferably operates at a sampling rate of at 
least twice the highest frequency of interest as determined by the digital filter 1620. In one 
embodiment, the sample rate compression block 1 622 reduces the sample rate by a factor of R, = 
37, corresponding to the number of samples during the period t as illustrated in Figure 10. The 



WO 98/46125 



PCT/US98/07131 



-27- 

output of the sample rate compression block 1622 provides one sample per time period t and 
thus four samples per time period T. The output of the sample rate compression block 1622 is a 
signal MF(k) (where k is a discrete index) which comprises approximately 1,266 samples per 
second. 

The signal MF(k) is provided as a first input to a first mixer 1624. The signal MF(k) is 
also provided as a first input to a second mixer 1626. A first demodulating signal £>,(*) is 
provided as a second input to the first mixer 1624, and a second demodulating signal D 2 (k) is 
provided as a second input to the second mixer 1626. The output of the first mixer 1624 is 
provided as an input to a first lowpass filter 1630, and the output of the second mixer is provided 
as an input to a second lowpass filter 1640. The bandwidths of the lowpass filters 1630, 1640 are 
preferably approximately 10 Hz. The signal MF(k) is also provided as a first input to a noise 
channel mixer 1628. A noise demodulating signal D 0 (k) is provided as a second input to the 
noise channel mixer 1628. The output of the noise channel mixer 1628 is provided to an input of 
a low pass filter 1650. The output of the low pass filter 1650 is provided to a sample rate 
compression block 1652. The output of the sample rate compression block 1652 is an estimate 
of the noise n(t). The output of the lowpass filters 1630 is provided to an input of a sample rate 
compressor 1632 and the output of the lowpass filter 1640 is provided to an input of a sample 
rate compressor 1642. The lowpass filter 1630 and the sample rate compressor 1632 together 
comprise a decimator 1631. The lowpass filter 1640 and the sample rate compressor 1642 
together comprise a decimator 1 64 1 , 

The output of the decimator 1631 is a signal S x (k) , which, as discussed below, is an 

estimate of the signal The output of the decimator 1641 is a signal S 2 (k) , which, as 

discussed below, is an estimate of the signal S 2 (k). As will be shown below, the selection of the 
first demodulating signal D\(k) and the second demodulating signal D 2 (k) in accordance with the 

present invention can reduce or eliminate the effects of noise in the two output signals 5, (it) 
and 5*2 (k) and also reduce or eliminate crosstalk between the two signals. 

The decimators 1632, 1642 decimate by a decimation rate R 2 . In a preferred 
embodiment, the decimators 1632, 1642 decimate by a decimation rate R2 = 20 to a sample rate 
of, for example, 63.3 Hz to provide a decimated output which can be further processed in 
accordance with the methods and apparatuses described in the above-referenced patents. The 
decimations which occur in the decimators 1632, 1642 reduce the rate at which the output signals 

S, {k) and S 2 (k) need to be processed while maintaining the sample rate well above the 10 Hz 
frequency content of the signals of interest. The outputs of the decimators 1632, 1642 are 
provided on respective output lines 1634 and 1644. 
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Decimating the signal MF(k) prior to demodulation, although not an approximation 
technique, can be simplified by assuming that each desired signal 5,(0 does not change 
appreciably during each period t. In many applications it is reasonable to assume that the desired 
signals 5,(0 and S 2 (t) will not change significantly during the time interval t shown in Figure 2. 

5 One skilled in the art will recognize that a sufficient condition for this assumption is that the 
highest significant frequency components in 5,(0 and $(/) are much lower than the modulation 
frequency. In the pulse-oximetry application the highest frequency of interest is typically around 
10 Hz, which is far below the 316.7 Hz fundamental of the modulation. Since n{t) is not a 
desired signal, no such assumption is necessary for «(/). Thus, while n(r) may vary erratically 

) over a modulation cycle, the signals 5,(0 and Sj(t) do not. Therefore, it is possible to perform 
pre-demodulation decimation that has little effect on 5,(0 and 5 2 (0 but may shape n(t) into n\t). 
The measured signal is decimated by a factor /?, = Q (where Q is the number of samples in a time 
period x) and then demodulated. 

Assuming R } = Q, then the spectral domain representation of the signal MF(k) at the 

i output of the sample rate compression block 1622 is given by (approximately): 

Since the sample rate compression block 1622 decimates at the same rate as the number 
of samples per period x, the decimation removes any x dependence in the expression for MF(f). 
The frequency components indexed by m increase four times faster than the frequency 
components indexed by n. This occurs because the modulated signals 5,(0 and 5 2 (0, which are 
indexed by n, occur in only one fourth of the samples, but the noise n(t), which is indexed by m, 
occurs in every sample. 

The demodulation operation can be performed either in the frequency or the time 
domain. A method for frequency domain demodulation of the signal MF(k) can be obtained by 
rewriting Equation 35 as: 



MF(j) =... MF_ 2 (j) + A/F, (/)+ MF 0 (f)+ MF x (f)+ MF 2 (ft+. . 



(36) 
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A/F 2 (/) = [5 I (/) + S 2 (/)]/r 

MF_ x {f) = {Sy(f)-S 2 {f)}IT 

MF Q (/) = [S { (/) + S 2 (/) + 4n'(f)} I T 

MF,V) = [S x {f)-S 2 {f)]IT (37) 

Aff r 2 (/) = [5 l (/)+5 2 (/)]/r 

A« r 3(/) = [5 l (/)-5 2 (/)]/r 

MF 4 (/) = [S x (/) + S 2 (/) + 4*'(/)] / r 



Where ri(k) is the decimated noise signal »(/). Estimates for the signal can be obtained by 
shifting the spectra of MF,(/ ) and MF 2 (/) by - 1/Fand -2/r, respectively, and then dividing the 
sum of the resultant by 2. Likewise, S 2 (f) can be obtained by dividing the difference of the 
5 resultant spectra by 2. In other words: 

S x ( / ) = MF X ( / - 1 / T) + MF 2 ( / - 2 / T) 
S 2 (f) = MF x {f-\IT)-MF 2 (f-2!T) 

Demodulation in the time domain is a more elegant method for obtaining S\(k) and S 2 (k). 
Time domain demodulation is obtained by using the frequency shift property of the Fourier 
transform given by: 

Ffco + oO <=> e iml f{t) (39) 

According to Equation 39, the frequency domain terms MFj(/) are related by a time shift in the 
10 time domain and this property can be used to generate the demodulation sequences D 0 -D 2 . A 
more complete development of this process (for the general case of N channels) is provided in 
Equations 42-50 below and in the text accompanying those equations. For the present case, 
where N= 2, using equations 42-50 gives: 

D 0 (k) = 0, I ft /, ... 
D } (k) = 1, -0.5, 0, -0.5, ... (40) 
D 2 (k) = 0, -0.J, 7, -0.5, ... 

The sequences shown in Equation 40 are repeating sequences of the four values shown. 
15 Thus, the demodulation waveforms are no more than short repeating sequences of simple 
coefficients. Since the samples MF(k) are time domain sequences, demodulation simply involves 
multiplying the samples MF(k) by the sequences in Equation 40. For example, the sequence of 
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coefficients A>(*) = (0, 1, 0, 1, ...) is provided to the multiplier 1628 to demodulate the signal 
MF{k) and produce the estimate of n(k). Similarly, the sequence of coefficients £»,(*) = (1, -0.5, 
0, -0.5, ...) is provided to the multiplier 1624 to demodulate the signal MF(k) and produce the 
estimate of £,(&). 

Multiple Channe l Modulation and Demodulation 

The two-channel pre-demodulation decimation technique described in the previous 
section can be extended to multi-channel systems having more than two desired signals. Figure 
17 illustrates an expansion of the two-channel modulator into a multi-channel 
modulator/demodulator. Figure 17 shows the first modulator 191 and the second modulator 193 
as shown in Figure 16. Further, Figure 17 shows a third modulator 1701 and an N*> modulator 
1 702. The signal input and a modulation input are provided to the first modulator 191 . 
The signal input S 2 (t) and a modulation input A/ 2 (0 are provided to the second modulator 193. A 
signal input 5 3 (/) and a modulation input M(0 are provided to the third modulator 1701. A 
signal input S„(t) and a modulation input Af*</) are provided to the modulator 1 702. 

The photodetector 150 is modeled as an adder 194 and an adder 196. The outputs of the 
modulators 191, 193, 1701, and 1703 are added together in the adder 194, to generate a 
composite signal M(t) where: 

M(t) = 5, (OA/, (0 + S 2 (t)M 2 (0 + S 3 (t)M J (t)+... +S N (t)M N (/) (41) 

The signal M{t) from the adder 194 is provided to the adder 196 where the signal M(t) is 
added to the signal n(t) which represents a composite noise signal caused by ambient light, 
electromagnetic pickup, and the like, which are also detected by the photodetector 150. The 
output of the adder 196 is the signal M(t) = A/(0 + «('), which includes the noise components as 
well as the signal components. 

The M(t) signal output of the adder 196 (i.e., the output of the detector 150) is applied to 
the input of the signal-processing block 1700. Within the signal-processing block 1700, the 
25 signal M(t) is first passed through an amplifier 197 and then through the analog bandpass filter 
198. The analog bandpass filter 198 provides anti-aliasing and removal of low frequency noise 
and DC. The desired signal components in the signals 5,(0 are frequency shifted by the operation 
of the modulation signals M{t) and are passed by the analog bandpass filter 198. 

The output of the analog bandpass filter 198 is sampled by the analog-to-digital 
converter 199 and converted therein to digital signals and provided to an input of the lowpass 
digital filter 1620. Output signals from the digital filter 1620 are provided to a sample rate 
compression block 1622, which reduces the sample rate by a decimation factor J?,. Together, the 
digital filter 1620 and the sample rate compression block 1622 comprise a decimator 1621. The 
output of the sample rate compression block 1622 is a signal MF{k). The signal MF{k) is 
provided as: the first input to the first mixer 1624; the first input to the second mixer 1626; a first 
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input to a third mixer 1710; a first input to an N th mixer 1712; and a first input to a noise 
channel mixer 1713. A first demodulating signal D\(k) is provided as a second input to the first 
mixer 1624. A second demodulating signal D 2 (k) is provided as a second input to the second 
mixer 1626. A third demodulating signal D£k) is provided to the third mixer 1710. A fourth 
5 demodulating signal D^k) is provided to the N* 1 mixer 1712. A noise demodulating signal D 0 (k) 
is provided to the noise channel mixer 1713. The outputs of the mixers 1624, 1626, 1710, 1712, 
and 1713 are provided as respective inputs of the lowpass filters 1630, 1640, 1720, 1730, and 
1740. The outputs of the lowpass filters 1630, 1640, 1720, 1730, and 1740 are provided as 
respective inputs of the decimators 1632, 1642, 1721, 1731 and 1741. Each of the decimators 
10 1632, 1642, 1 721, 173 1 and 1741 reduces the sample rate by a decimation rate R 2 . 

The output of the sample rate compressor 1632 is a signal S l (k) , which, as discussed 
below, is an estimate of the signal S\(k). Likewise, the output of the sample rate compressor 
1642 is an estimate of S 2 (t) t the output of the sample rate compressor 1721 is an estimate of the 
signal 5 3 (/), the output of the sample rate compressor 1731 is an estimate of the signal S u (t\ and 
1 5 the output of the sample rate compressor 1 741 is an estimate of the signal n{t). 

As will be shown below, the selection of the demodulating signals D,(i) for / = 0..JV in 
accordance with the present invention can substantially reduce or eliminate the effects of noise in 

the output signals S f (k) and n(k\ and can also substantially reduce or eliminate crosstalk 
between the signals. 

20 As shown in Figure 1 7, a set of N+l signals S f [k] i = 1 JV, and n(k) are sampled at a rate 

T/QN, where T is a modulation period. For simplicity, the decimation rate R^ is assumed to be 
the same as the factor Q. The assumption that R\ = Q is not a necessary assumption, but rather is 
used here to simplify the mathematics. The signals are combined according to the formula: 



S(k)= Mi(k)Si(k)+ M 2 (k)S 2 (k)+ Ms(k)Ss(k)+ ... + M N (k)S N (k) + n(k) 
Using the symbol * to denote the convolution operator, the terms M,(k) are given by: 



(42) 



25 where 



HW = A(2M/D*/>,(0 \ imkriQN 

M t (k) = A(2M/D */>(/) \ imkT/QN (43) 
AM*) = A(2M/ T)*P N (0 \,. kTIQN 



A ( X )=[ 1 ,:c| - a5 44 

|0 otherwise 
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and 



P,(0= £ h(t-nT) 



(45) 



(where 8(A) is the Kroneker delta function, which is 1 for* « 0, and 0 for all other values of*), 
and 



N 



(46) 



After the pre-demodulation and sample rate compression stage 1622, which decimates by a factor 
'• Q, the signal in the frequency domain is given approximately by 



where 



n ^ j2nkn 

s* ~ e n 



The demodulator sequences are then given by: 



2 



(47) 



(48) 



£,(*) = />(*)- 



Do(k) 
N 



N 



(49) 
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where 

Pi(k) = Pi(t)[^ t p 2 (k) = P 2 (t)\ l _*L,,.. t P N (k)=P N (t)\ im kT (50) 

2" 2N IN 

The post demodulation lowpass filters 1630, 1640, 1720, 1730 and 1740, and the post 
demodulation sample rate compression stages 1632, 1642, 1721, 1731 and 1741 suppress high 
frequency artifacts which are produced by the modulation/demodulation process. Note that 
5 Equation 49 reduces to Equation 40 for N = 2. 
Adaptive Demodulation 

The multi-channel pre-demodulation decimation technique described in the previous 
section can be extended to an adaptive multi-channel system having an adjustable pre- 
demoduiation decimation rate and an adjustable post-demodulation decimation rate. Figure 18 
10 illustrates an expansion of the multi-channel modulator into a adaptive multi-channel 
modulator/demodulator 1800. Figure 18 shows the first modulator 191 and the N 01 modulator 
1702 as shown in Figure 17. The signal input S x (t) and a modulation input M,(r) are provided to 
the first modulator 19 1 . A signal input 5 N (/) and a modulation input A/ N (0 are provided to the N* 1 
modulator 1702. 

15 The photodetector 150 is modeled as an adder 194 and an adder 196. The outputs of the 

modulators 191, 193, 1701, and 1703 are added together in the adder 194, to generate a 
composite signal M(t) where: 

M(t) = S x (t)M t (0 + . + S N (t)M N (0 (51) 

The signal M(t) from the adder 194 is provided to the adder 196 where the signal M(t) is 
added to the signal n(t) which represents a composite noise signal caused by ambient light, 

20 electromagnetic pickup, and the like, which are also detected by the photodetector 150. The 
output of the adder 196 is the signal M(t) = M(t) + n(t), which includes noise components as well 
as the signal components. 

The M(t) signal output of the adder 196 (i.e., the output of the detector 150) is applied to 
the input of the signal-processing block 1800. Within the signal-processing block 1800, the 

25 signal Af(/) is first passed through the amplifier 197 and then through the analog bandpass filter 
198. The analog bandpass filter 198 provides anti-aliasing and removal of low frequency noise 
and DC. The desired signal components in the signals S t (t) are frequency shifted by the operation 
of the modulation signals M,{t) and are passed by the analog bandpass filter 198. 

The output of the analog bandpass filter 198 is sampled by the analog-to-digital 

30 converter 199 and converted therein to digital signals and provided to an input of a decimation 
block 1820. The adaptive decimation block 1820 comprises a digital lowpass filter and a sample 
rate compressor that reduces the sample rate by the decimation rate /?,. The filter coefficients 



WO 98/46125 



PCT/US98/07131 



10 



15 



-34- 

and decimation rate rt, are provided to a control input of the adaptive decimation block 1820 by 
an output of an adaptive algorithm block 1850. Equation 35 assumes that the decimation rate R { 
is equal to Q, However, in general, the value of Q may be different than the decimation rate J?,. 
The output of the adaptive decimation block 1820 is a signal MF(k). 

The signal MF{k) is provided to the first input of the first mixer 1624, to the first input of 
the N 'th mixer 1712, and to the first input of the noise channel mixer 1713. A first demodulating 
signal D\(k) is provided to a second input of the first mixer 1624 from a signal generator 1841. 
The fourth demodulating signal D^k) is provided to the mixer 1712 from an output of a 
signal generator 1831. The noise demodulating signal D 0 (k) is provided to the noise channel 
mixer 1713 from an output of a signal generator 1832. A control input to each of the signal 
generators 1831, 1832, and 1841 is provided by the output of the adaptive algorithm 1850. In yet 
another embodiment, the adaptive algorithm 1850 may also be controlled by other signal 
processing elements downstream of the signal processor 1800. 

The outputs of the mixers 1713, 1624, and 1712 are provided as respective inputs to 
adaptive decimation blocks 1840, 1830, and 1834 respectively. Each of the adaptive decimation 
blocks 1 840, 1 830, and 1 834 has a control input provided by the output of the adaptive algorithm 
block 1 850. The output of the adaptive decimation block 1 840 is an estimate of the signal n{t) 
and it is provided to an input of the adaptive algorithm block 1850. In an alternate embodiment, 
the signal estimates S f (k) are also provided to the adaptive algorithm block 1 850. 

An output of the decimator 1830 is a signal S t (k) , which, as discussed above, is an 
estimate of the signal £,(*). Likewise, the output of the decimation block 1834 is an estimate of 
the signal SJf). As shown above, the selection of the demodulating signals D,(i) for i = 0..JV in 
accordance with the present invention substantially reduces or eliminates the effects of noise in 

the output signals S^k) and n(k\ and also substantially reduces or eliminates crosstalk between 
25 the signals. 

As shown in Figure 18, a set of N+l signals Sflt] i = 1 ..JV, and n{k) are sampled at a rate 
TIQN, where T is a modulation period, and R x is the decimation rate of the decimation block 
1 820. The signals are combined according to the formula: 

S(k) = Mi(k)Si(k) + . . . + M N (k)S N (k) + n(k) (5 2) 



20 



30 



Each of the adaptive decimators 1820, 1840, 1830, and 1 834 comprises a digital lowpass 
filter and a sample rate compressor. The characteristics of the digital lowpass filters (e.g., the 
number of filter coefficients and values of the filter coefficients) and the sample rate compression 
factor of each adaptive decimator is provided to a control input of the adaptive decimator. The 
control inputs are driven by an adaptive algorithm 1850. The signal generators 1831, 1832 and 
1841 generate the demodulation sequences for the demodulators 1624, 1712, and 1713 
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respectively. The demodulation sequences produced by the signal generators 1831, 1832 and 
1 84 1 are controlled by the adaptive algorithm 1 850. 

The adaptive algorithm adjusts the pre-demodulation decimation rate R x (in the adaptive 
demodulator 1820), and the post-demodulation decimation rate R 2 (in the adaptive demodulators 
5 1830, 1834 and 1840) according to the noise in the noise estimate n(k) 1746 and (optionally) 
according to the signals S f (k) . The product R X R 2 is the total decimation rate from the signal 

S(k) at the output of the A/D converter 199 to the signals S,(£)at the output of the signal 
processing block 1800. The adaptive algorithm may adjust R x and R 2 such that the product R } R 2 
varies, or the adaptive algorithm may adjust /?, and R 2 such that the product R X R 2 is substantially 
1 0 constant Typically, the adaptive algorithm will keep the R X R 2 product constant so that the signal 
processing blocks downstream of the signal processor 1800 will operate at a substantially 
constant sample rate. 

Typically, each of the signal generators 1841, 1831 and 1832 generates a repeating 

sequence of numbers. The number of elements in the sequence is a function of the decimation 
15 factor R { . As discussed above in connection with Figure 3, when /?, = I, there are preferably 

148 values in each demodulation sequence. As discussed above in connection with Figure 17, 

when R\ = 37, there are preferably only 4 values in the demodulation sequences. 

The adaptive algorithm selects R u R 2 , and the filter transfer functions in the adaptive 

decimators 1820, 1830, 1834, and 1840 to improve the quality of the output signals S f (k) . For 
20 example, in high ambient noise environments, the higher order harmonics of the output signals 

are often contaminated by ambient noise (as discussed in connection with Figures 14 and 20). 

Thus, the higher order harmonics are preferably not demodulated when ambient noise is present. 

To avoid demodulation of the higher order harmonics the adaptive demodulator 1850 can set 

= I and R 2 = 37, and thereby demodulate according to the method described in connection with 
25 Figures 3-14. Alternatively, the adaptive demodulator 1850 can set R x = 37, set R 2 = 1, and set 

the transfer function of the lowpass filter in the adaptive decimator 1 820 to provide a very fast 

rolloff (thereby filtering out the higher order harmonics). 

Conversely, in low ambient noise environments, the higher order harmonics of the 

output signal are less contaminated by ambient noise, and thus the higher order harmonics may 
30 be demodulated. In one embodiment, to demodulate the higher order harmonics, the adaptive 

demodulator 1850 can set R } = 37 and set R 2 « 1, to demodulate according to the method 

described in connection with Figure 17. This is especially advantageous when perfusion is low, 

because, when perfusion is low the output signals S f (k) are typically very weak and are 
contaminated by random noise. Demodulating more of the higher order harmonics increases the 
35 signal-to-noise ratio because it adds the harmonics (which are correlated) to the output signals, 
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and tends to average out the noise (which is uncorreiated). Thus, the signal strength increases, 
and the noise is reduced. 

One skilled in the art will recognize that the examples in the preceding two paragraphs 
are merely two points on a continuum and that the adaptive algorithm 1850 can generate many 
desirable solutions on the continuum. 
Ambient Light Rejection 

In the pulse oximeter, one of the major contributors to the noise signal «(/) is ambient 
light that is detected by the photodetector 150. One aspect of the present invention 
advantageously provides a method for choosing the modulation sampling rate / s and the factor Q 
so that the effects of ambient light can be removed by the post demodulation filtering and 
decimation stages. Note that Q is the number of samples during the on period (i.e., modulation 
signal sample turn on time Q) and is preferably also the decimation rate R x for the pre- 
demodulation sample rate compressor 1 622 (in general the values of Q and R x may be different). 
The particular embodiment described by Equation 35 assumes that the value Q is also used as 
decimation rate for the pre-demodulation decimator 1820. 

In the system shown in Figures 3 and 16, which demodulates two harmonics, the period 
of a modulation cycle is given by: 



T = 4Q/f, 

where/ s is the sample rate. Defining the two line equations 



where 



(53) 



(54) 



f a = line frequencies of concern 
n = line frequency harmonic numbers of concern ^ 



then the effects due to ambient light will be minimized when 



\y(f a , n)\>SBF 

\z{f a ,ri)\>SBF (56) 



where SBF is the stop band frequency of the post demodulation and decimation stages (e.g., the 
10 Hz iowpass filter 1630 and the sample rate compressor 1632, etc.). 
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Figure 19 is a flowchart showing a method for selecting^ and O. The method begins at 
a process block 1902 wherein the ambient light frequencies / a and important harmonic 
components n are identified. Important harmonics are defined as those harmonics that will 
degrade system performance below acceptable levels when detected by the detector 150. The 
5 process then advances from the process block 1902 to a process block 1904. In the process block 
1904, the values of/ a and n identified in the process block 1902 are used in conjunction with 
Equation 54 to identify a collection of acceptable values of T. Upon completion of the process 
block 1904, the process advances to a process block 1806. In the process 1906, suitable values of 
/ s and Q are chosen using the values of IT obtained in the process block 1904 and the equation T= 
10 *Q!U One skilled in the art will recognize that, since T is proportional to the ratio of Q/f S9 
knowing Twill not uniquely determine either / s or Q. 

For example, given power line frequencies of 50 ± I Hz and 60 ± 1 Hz then the range of 
/ a is given by approximately the union of the interval 49-51 Hz and the interval 59-61 Hz, which 
can be expressed mathematically as: 

f a <*[49 t 51]v[59,61] (57) 

15 Assuming that all harmonics up to the 18 th harmonic are to be suppressed, then n = 1 ...18. In a 
preferred embodiment, using these values for / a and «, application of the method in Figure 19 
results in/ s = 46,875 Hz and acceptable Q values of 37 and 4 1 . 

The process leading to Equation 57 is illustrated graphically by Figure 20, where the 
harmonics of the ambient light frequency^ (in Hz) are plotted versus the plethysmograph 

20 signal frequency (also in Hz). Figure 20 has an x-axis showing the ambient light frequency 
from 44 Hz to 64 Hz. The ambient light frequency will usually correspond to the frequency 
of the power lines, which is nominally 60 Hz (in the U.S.) and 50 Hz (outside the U.S.). 
However, power line frequency regulation typically varies somewhat, and thus Figure 20 
shows frequencies above and below the nominal frequencies. 

25 Figure 20 also shows a y-axis showing the plethysmograph signal frequency from -10 

Hz to 10 Hz. One skilled in the art will recognize that negative frequencies occur in the 
mathematics described above. In particular, a signal that is modulated from baseband up to some 
carrier frequency will exhibit two sidebands, a sideband above the carrier frequency 
corresponding to the frequency of the baseband signal, and a sideband below the carrier 

30 frequency corresponding to the negative of the baseband frequency. Thus, when dealing with 
modulation and demodulation, it is convenient to deal with positive and negative frequencies. 

Figure 20 also shows harmonic lines corresponding to the 5 th , 6 th , 7 th , 10 th , 1 1 th , 12 th , 13 th , 
and 14 th harmonics of the ambient light frequency. The harmonic lines correspond to the 
harmonics produced in the plethysmograph signal by the demodulation (mixing down) of 

35 harmonics of the power line frequency. The lines in Figure 20 are calculated using Equation 54 
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for 1/r = 3 16.72 Hz. Some of the harmonic lines correspond to y{f w n) and some correspond to 
z(fc»n) from Equation 54. Harmonic lines that are not shown (e.g., the line corresponding to the 
8 th harmonic) fall outside the displayed limits of the x-axis and y-axis. 

Figure 20 can be used to determine the stop band frequencies as shown in Equation 56. 

5 For example, the harmonic lines in Figure 20 show that for an ambient light frequency of 49 Hz, 
the 13 th harmonic of the ambient light frequency will appear in the plethysmograph signal at 
approximately 3 Hz. Thus, Figure 20 shows that for plethysmograph bandwidth of 10 Hz, none 
of the first 14 harmonics of the ambient light will appear in the plethysmograph signal for 
ambient light frequencies between approximately 61.2 Hz and approximately 58.5 Hz, which is 

10 consistent with Equation 57. The first ambient harmonics that do appear for a plethysmograph 
bandwidth of 10 Hz are the 5 th harmonic and the 11 * harmonic. 
Other Embodiments 

In the preferred embodiment of the present invention, the hardware described above 
is implemented in a digital signal processor and associated circuitry. The LED modulation 

15 block 104 and the LED demodulation state table block 352 comprise algorithms implemented 
by program code executed by the digital signal processor. In addition, the configuration 
variables, such as for example, the hardware delay value, the hardware distortion value and 
the hardware scale value are provided as inputs to the digital signal processor when it is set 
up. For example, the main operating program of the digital signal processor may be stored in 

20 non-volatile ROM or PROM, and the variables may be stored in flash memory during a setup 
procedure. Techniques for communicating to and from a digital signal processor during such 
setup procedures are well known to persons of skill in the art, and will not be described in 
detail herein. For example, the configuration bus 310, discussed above, represents a 
communication path to the flash memory during such a setup procedure. The data provided to 

25 the configuration bus 3 10 may be provided by a system operator (not shown) or the data may 
be provided from look-up tables (not shown) maintained for different embodiments of the 
LEDs 106, 108 and the detector 150. 

Although described above in connection with a pulse oximetry system wherein a 
parameter to be measured is the attenuation of red and infrared light passing through a portion of 

30 a subject's body, it should be understood that the method and apparatus described herein can also 
be used for other measurements where two or more signals are passed through a system to be 
analyzed. In particular, the present invention can be used to demodulate two combined 
parametric signals responsive to the system to be analyzed where the two parametric signals have 
a predetermined timing relationship between them, as described herein. 

35 One skilled in the art will recognize that the lowpass filters provided in connection with 

the decimation blocks may provide other filter functions in addition to lowpass filtering. Thus, 
for example, the lowpass filters 1620, 1622, 1630, 1640, 1650, 1720, 1730, and 1740, and the 
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decimators 1820, 1830, 1834, and 1840 may provide other filter functions (in addition to lowpass 
filtering) such as, for example, bandpass filtering, bandstop filtering, etc. Moreover, the post- 
demodulation decimation rate R 2 need not be the same for each output channel. Thus, for 
example, in Figure 18, the decimator 1840 may have a first decimation rate R 2 = r, while the 

5 decimators 1 830 and 1 834 have a second decimation rate R 2 = r 2 . 

Although described above in connection with a particular embodiment of the present 
invention, it should be understood the description of the embodiment is illustrative of the 
invention and are not intended to be limiting. Various modifications and applications may occur 
to those skilled in the art without departing from the true spirit and scope of the invention as 

10 defined in the appended claims. 
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WHAT IS CLA1MHD IS : 



1 1. 



2 



11 
12 
13 
14 
15 



An apparatus for measuring blood oxygenation in a subject, said apparatus comprising: 
a first signal source which applies a first input signal during a first time interval; 

3 a second signal source which applies a second input signal during a second time 

4 interval; 

5 a detector which detects a first parametric signal responsive to said first input 
signal passing through a portion of said subject having blood therein and which detects a 
second parametric signal responsive to said second input signal passing through said 
portion of said subject, said detector generating a detector output signal responsive to 

9 said first and second parametric signals; and 

10 a signal processor which receives said detector output signal, said signal 

processor demodulating said detector output signal by applying a first demodulation 
signal to a signal responsive to said detector output signal to generate a first output signal 
responsive to said first parametric signal and applying a second demodulation signal to 
said signal responsive to said detector output signal to generate a second output signal 
responsive to said second parametric signal, each of said first demodulation signal and 

16 said second demodulation signal comprising at least a first component having a first 

17 frequency and a first amplitude and a second component having a second frequency and 

1 8 a second amplitude, said second frequency being a harmonic of said first frequency, said 

19 second amplitude selected to be related to said first amplitude to minimize crosstalk 

20 fr° m said first parametric signal to said second output signal and to minimize crosstalk 

21 from said second parametric signal to said first output signal. 

2. The method of Claim 1 , wherein the second amplitude is determined by turning 
2 off one of the first and second signal sources and measuring the crosstalk between one of the 

parametric signals and the non-corresponding output signal while varying the second amplitude 

4 and selecting a second amplitude which minimizes the measured crosstalk. 

3. A method of minimizing crosstalk between two signals generated by applying a 
2 first pulse and a second pulse to measure a parameter, wherein said first pulse and said second 

pulse are applied periodically at a first repetition rate defining a period, and wherein said first 
pulse is generated during a first interval in each period and said second pulse is generated during 

5 a second interval in each period, said second interval spaced apart from said first interval, said 
first and second pulses producing first and second parametric signals responsive to said 
parameter, said first and second parametric signals being received by a single detector which 
outputs a composite signal responsive to said first and second parametric signals, said method 

9 comprising the steps of: 
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10 applying a first demodulation signal to said composite signal to generate a first 

11 demodulated output signal, said first demodulation signal comprising at least a first 

12 component having a first frequency corresponding to said first repetition rate and having 

13 a first amplitude, said first demodulation signal further comprising a second component 
1^ having a second frequency which is a harmonic of said first frequency and having a 

15 second amplitude which has a selected proportional relationship to said first amplitude; 

16 applying a second demodulation signal to said composite signal to generate a 

17 second demodulated output signal, said second demodulation signal comprising said first 

18 component at said first frequency and said first amplitude and comprising said second 

19 component at said second frequency and said second amplitude, at least one of said first 

20 and second components of said second demodulation signal having a selected phase 

21 difference with respect to the corresponding one of said first and second components of 

22 said first demodulation signal; and 

23 lowpass filtering said first demodulated output signal to generate a first 

24 recovered output signal responsive to said first parametric signal; and 

25 lowpass filtering said second demodulated output signal to generate a second 
26. recovered output signal responsive to said second parametric signal. 

1 4. The method as defined in Claim 3, wherein said selected phase difference is 7t. 

1 5. The method as defined in Claim 3, wherein: 

2 said first pulse and said second pulse are generally rectangular pulses having a 

3 duty cycle, and wherein said rectangular pulses comprise a plurality of sinusoidal 

4 components including a fundamental component corresponding to said first frequency 

5 and a first harmonic component corresponding to said second frequency, said 

6 fundamental component having a fundamental component amplitude and said first 

7 harmonic component having a first harmonic component amplitude, said first harmonic 

8 component amplitude related to said harmonic component amplitude by a first 

9 proportionality value; and 

1 0 said second amplitude of said second component of said first demodulation 

1 1 signal is related to said first amplitude of said first component of said first demodulation 

12 signal by a second proportionality value which is approximately the inverse of said first 

1 3 proportionality value. 

1 6. The method as defined in Claim 3, further including the steps of: 

2 sampling said composite signal when neither said first pulse nor said second 

3 pulse is active to obtain a sampled signal; and 

4 measuring said sampled signal to determine a noise level of said parametric 

5 signals. 

1 7. The method as defined in Claim 3, further including the steps of: 
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performing a transform on said composite signal to generate a spectra of said 
composite signal; 

sampling said spectra at a plurality of frequencies other than at predetermined 
ranges of frequencies around said first frequency and around harmonics of said first 
frequency; 

determining an average of the magnitudes of said sampled plurality of 
frequencies; and 

comparing said average to a selected threshold to determine whether the average 
magnitude exceeds said selected threshold. 

8. A method of demodulating a composite signal generated by applying first and 
second periodic pulses of electromagnetic energy to a system having a parameter to be measured 
and by receiving signals responsive to said electromagnetic energy after having passed through 
said system and being affected by said parameter being measured, said signals received as a 
composite signal having components responsive to said first and second pulses, said method 
comprising the steps of: 

applying a first demodulation signal to said composite signal to generate a first 
demodulated signal, said first demodulation signal comprising a first component having 
a first frequency corresponding to a repetition frequency of said first and second pulses 
and comprising a second component having a frequency which is a harmonic of said first 
frequency, said first component having a first amplitude and said second component 
having a second amplitude, said second amplitude having a predetermined relationship 
to said first amplitude, said predetermined relationship selected to cause said first 
demodulated signal to have low frequency components responsive only to said first 
pulse; and 

lowpass filtering said first demodulated signal to generate a first output signal, 
said first output signal varying in response to an effect of said parameter on the 
electromagnetic energy received from said first pulse. 

9. The method as defined in Claim 8, further including the steps of: 

applying a second demodulation signal to said composite signal to generate a 
second demodulated signal, said second demodulation signal having first and second 
components corresponding to said first and second components of said first 
demodulation signal, at least one of said first and second components of said second 
demodulation signal having a selected phase relationship with the corresponding one of 
said first and second components of said first demodulation signal; and 

lowpass filtering said second demodulated signal to generate a second output 
signal, said second output signal varying in response to an effect of said parameter on the 
electromagnetic energy received from said second pulse. 
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1 1 0. The method as defined in Claim 9, wherein said selected phase relationship is a 

2 n phase difference. 

1 1 1 . A pulse oximetry system, comprising: 

2 a modulation signal generator, said modulation signal generator generating a 

3 first modulation signal comprising a first pulse which repeats at a first repetition 

4 frequency, said first pulse having a duty cycle of less than 50%, said modulation signal 

5 generator generating a second modulation signal comprising a second pulse which also 

6 repeats at said first repetition frequency, said second pulse having a duty cycle of less 

7 than 50%, said second pulse occurring at non-overlapping times with respect to said first 

8 pulse, said first and second pulses comprising a plurality of components wherein a first 

9 component has a frequency corresponding to said repetition frequency and a second 

10 component has a second frequency corresponding to twice said first frequency, said 

11 second component having an amplitude which has a first predetermined relationship to 

1 2 an amplitude of said first component; 

13 a first transmitter which emits electromagnetic energy at a first wavelength in 

14 - response to said first pulse; 

15 - a second transmitter which emits electromagnetic energy at a second 

1 6 wavelength in response to said second pulse; 

17 a detector which receives electromagnetic energy at said first and second 
19 <■ wavelengths after passing through a portion of a subject and which generates a detector 

19 output signal responsive to the received electromagnetic energy, said detector output 

20 signal including a signal component responsive to attenuation of said electromagnetic 

21 energy at said first wavelength and a signal component responsive to attenuation of said 

22 electromagnetic energy at said second wavelength; 

23 a first demodulator which multiplies said detector signal by a first demodulation 

24 signal and generates a first demodulated output signal, said first demodulation signal 

25 comprising a first component having said first frequency and having a first amplitude 

26 ar >d comprising a second component having said second frequency and having a second 

27 amplitude, said second amplitude having a second predetermined relationship to said 

28 first amplitude which second predetermined relationship is inversely proportional to said 

29 first predetermined relationship; and 

30 a second demodulator which multiplies said detector signal by a second 

31 demodulation signal and generates a second demodulated output signal, said second 

32 demodulation signal comprising a first component having said first frequency and 

33 having said first amplitude and comprising a second component having said second 

34 frequency and having said second amplitude, at least one component of said second 
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35 demodulation signal having a selected phase relationship with a corresponding one 

36 component of said first demodulation signal. 

1 12. The method as defined in Claim 11, wherein said selected phase relationship is a 

2 n phase difference. 

1 13, A method of minimizing crosstalk between two signals generated by applying a 

2 first pulse and a second pulse to measure a parameter, wherein said first pulse and said second 

3 pulse are applied periodically at a first repetition rate defining a period, and wherein said first 

4 pulse is generated during a first interval in each period and said second pulse is generated during 

5 a second interval in each period, said second interval spaced apart from said first interval, said 

6 first and second pulses producing first and second parametric signals responsive to said 

7 parameter, said first and second parametric signals being received by a single detector which 

8 outputs a composite signal responsive to said first and second parametric signals, said method 

9 comprising the steps of: 

10 providing said composite signal to an analog to digital converter to produce a 

1 1 sequence of digital values; 

12 decimating said sequence of digital values to produce a decimated sequence of 

13 digital values; 

14 applying a first sequence of demodulation coefficients to said decimated 

15 sequence of digital values to generate a first demodulated output signal; 

16 applying a second sequence of demodulation coefficients to said decimated 

1 7 sequence of digital values to generate a second demodulated output signal; 

18 lowpass filtering said first demodulated output signal to generate a first 

19 recovered output signal responsive to said first parametric signal; and 

20 lowpass filtering said second demodulated output signal to generate a second 

21 recovered output signal responsive to said second parametric signal. 

1 14. The method as defined in Claim 13, further comprising the steps of: 

2 applying a third sequence of demodulation coefficients to said decimated 

3 sequence of digital values to generate a third demodulated output signal; and 

4 lowpass filtering said third demodulated output signal to generate a recovered 

5 output signal responsive to noise detected by said detector. 

1 1 5. A method of demodulating a composite signal generated by applying a plurality 

2 of periodic pulse trains of electromagnetic energy to a system having a parameter to be measured 

3 and by receiving signals responsive to said electromagnetic energy after having passed through 

4 said system and being affected by said parameter being measured, said signals received as a 

5 composite signal having components responsive to said plurality of periodic pulse trains, each of 

6 said components corresponding to one each of said pulse trains, said method comprising the steps 

7 of: 



WO 98/46125 



PCT/US98/07131 



-45- 

8 sampling said composite signal using an analog to digital converter to produce a 

9 sequence of digital values; 

10 decimating said sequence of digital values to produce a decimated sequence of 

1 1 values; 

12 providing said decimated sequence of values to a first input of each of plurality 

13 of multipliers; 

14 providing a plurality of sequences of demodulation coefficients to a second 

15 input of each of said multipliers such that each multiplier is provided a unique sequence 

16 of demodulation coefficients; and 

17 lowpass filtering the output of each of said multipliers such that the output 

18 signal of each of said lowpass filters approximately corresponds to one each of said 

19 components. 

1 1 6. A pulse oximetry system, comprising: 

2 a modulation signal generator, said modulation signal generator generating a 

3 first modulation signal comprising a first pulse which repeats at a first repetition 

4 frequency, said first pulse having a duty cycle of less than 50%, said modulation signal 

5 generator generating a second modulation signal comprising a second pulse which also 

6 repeats at said first repetition frequency, said second pulse having said duty cycle of less 

7 than 50%, said second pulse occurring at non-overlapping times with respect to said first 
1 8 pulse; 

9 a first transmitter which emits electromagnetic energy at a first wavelength in 

I o response to said first pulse; 

II a second transmitter which emits electromagnetic energy at a second 

1 2 wavelength in response to said second pulse; 

13 a detector which receives electromagnetic energy at said first and second 

14 wavelengths after passing through a portion of a subject and which generates a detector 

15 output signal responsive to the received electromagnetic energy, said detector output 

16 signal including a signal component responsive to attenuation of said electromagnetic 

17 energy at said first wavelength and a signal component responsive to attenuation of said 

1 8 electromagnetic energy at said second wavelength; 

19 a sampling analog to digital converter which converts said detector output into a 

20 sequence of digital values; 

21 a decimator which decimates said sequence of digital values to produce a 

22 decimated sequence; 

23 a first demodulator which multiplies said decimated sequence by a 

24 demodulation sequence and generates a first demodulated output signal; and 
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a second demodulator which multiplies said decimated sequence by a second 
demodulation sequence and generates a second demodulated output signal. 

17. The apparatus of claim 16 further comprising a third demodulator which 
multiplies said decimated sequence by a third demodulation sequence and generates a third 
demodulated output signal corresponding to noise produced by said detector. 

1 8. The apparatus of claim 1 6 wherein said decimator has a decimation rate equal to 
the number of said samples produced during said duty cycle. 

19. In a system that includes: a sampling frequency generator generating a sampling 
frequency / s ; a detector; a modulation signal generator, said modulation signal generator 
generating a sequence of pulses having a pulse repetition frequency, said pulses having a duty 
cycle of Q sample periods of said sampling frequency / s ; a transmitter which emits 
electromagnetic energy at a wavelength in response to said pulse, said detector receiving said 
electromagnetic energy to generate a detector output signal responsive to the received 
electromagnetic energy, said detector output signal including noise caused by ambient 
electromagnetic energy detected by said detector; a digital to analog converter which produces 
digital samples of said detector output signal at said sampling frequency; and a demodulator 
which demodulates said digital samples to produce a desired output signal, a method for 
minimizing said noise due to ambient electromagnetic energy in said desired output signal 
comprising the steps of: 

identifying all undesired frequency components of said ambient electromagnetic 
energy detected by said detector; 

using said undesired frequency components to compute a set of acceptable 
modulation cycle times T\ and 

using said acceptable modulation cycle times to select said / s and said Q using 
an equation T=4Q/f s . 

20. A method for demodulating a multi-channel composite signal comprising the 

acts of: 

generating a demodulation sequence for a selected channel; 
providing said demodulation sequence to a first input of a demodulator; 
providing a sampled composite signal to a second input of said demodulator. 

21 . The method of Claim 20, further comprising the act of decimating an output 
of said demodulator by a decimation factor R. 

22. The method of Claim 2 1 , wherein said act of decimating includes the acts of: 
lowpass filtering; and 

sample rate compressing. 
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23. The method of Claim 20, further comprising the act of decimating said sampled 
composite signal before providing said sampled composite signal to said second input of said 
demodulator. 

24. The method of Claim 23, wherein said act of decimating includes the acts of 
iowpass filtering and sample rate compressing. 

25. The method of Claim 24, wherein said acts of Iowpass filtering and sample rate 
compressing are controlled by an adaptive algorithm. 

26. The method of Claim 20, further comprising the acts of : 

decimating in a first decimator said sampled composite signal before providing 
said sampled composite signal to said second input of said demodulator, wherein said 
first decimator has a first Iowpass filter transfer function and a first decimation rate; and 

decimating in a second decimator an output of said demodulator, wherein 
said second decimator has a second Iowpass filter transfer function and a second 
decimation rate. 

27. The method of Claim 26, further comprising the act of controlling said first 
Iowpass filter, said second Iowpass filter, said first decimation rate, and said second decimation 
rate by an adaptive algorithm. 

28. The method of Claim 26, wherein a product of said first decimation rate and 
said second decimation rate is substantially constant. 
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